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Résumé 


Dans cette thèse, notre objectif principal est de développer un cadre analytique et numé- 
rique visant à définir le spectre quadratique généralisé des opérateurs bornés, qui est déjà 
défini pour les matrices. De plus, nous cherchons à étendre les recherches existantes sur le 
spectre généralisé de deux opérateurs bornés et son application dans l’approximation du 
spectre d'opérateurs non bornés. Cette généralisation est essentielle pour la création de la 
méthode du spectre quadratique généralisé, visant à résoudre le problème de la pollution 
spectrale, un phénomène couramment observé dans l’approximation spectrale des opérateurs 
non bornés. Ces opérateurs sont utilisés pour modéliser mathématiquement de nombreux 
problèmes dans divers domaines tels que la physique, la chimie, la mécanique quantique, 
l'optique, le traitement du signal, l’électromagnétisme, l’informatique et l’informatique quan- 


tique. 
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Classification Mathématique des Sujets (2010) : 35P15, 47A75, 35J10. 


Abstract 


In this thesis, our main objective is to develop an analytical and numerical framework 
to define the generalized quadratic spectrum of bounded operators, à concept already es- 
tablished for matrices. Furthermore, we aim to extend existing research on the generalized 
spectrum of two bounded operators and its application in approximating the spectrum of 
unbounded operators. This generalization is essential for the creation of the generalized qua- 
dratic spectrum method, designed to address the issue of spectral pollution, a phenomenon 
commonly observed in the spectral approximation of unbounded operators. These operators 
are used to mathematically model various problems in diverse fields such as physics, chemis- 
try, quantum mechanics, optics, signal processing, electromagnetism, computer science, and 


quantum computing. 
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Introduction Générale 
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1.1 Historique de la théorie spectrale 


L'histoire de la théorie spectrale remonte aux 18e et 19e siècles, et elle a évolué de manière 
complexe au fil du temps. Au 18e siècle, des mathématiciens comme Leonhard Euler et 
Joseph Louis Lagrange ont commencé à explorer les valeurs propres d'opérateurs linéaires, 
en se concentrant principalement sur des opérateurs différentiels. Au début du 19e siècle, 
Joseph Fourier à apporté d'importantes contributions en développant la série de Fourier 
pour analyser les signaux périodiques. La théorie spectrale à été appliquée à la résolution 
d'équations différentielles par des mathématiciens comme Augustin-Louis Cauchy et Jean- 
Baptiste Joseph Fourier, qui ont introduit des opérateurs différentiels pour étudier les modes 
de vibration des cordes et des plaques. 

Au milieu du 19e siècle, le mathématicien allemand Peter Gustav Lejeune Dirichlet a joué 
un rôle majeur en étudiant les problèmes aux valeurs propres pour les équations aux dérivées 
partielles, formulé des conditions aux limites portant son nom (conditions de Dirichlet). Au 
tournant du 20e siècle, David Hilbert a développé les espaces de Hilbert, qui sont essentiels 
en théorie spectrale et en mécanique quantique. 

Cette dernière a émergé au début du 20e siècle, avec des mathématiciens tels qu'Erwin 
Schrüdinger et Werner Heisenberg utilisant la théorie spectrale pour formuler les équations de 
la mécanique quantique et comprendre les propriétés spectrales des opérateurs d’observable. 

La théorie spectrale s’est continuellement développée au 20e siècle, avec des contributions 
importantes de mathématiciens comme John von Neumann. Elle a donné naissance à des 
concepts tels que le spectre discret, le spectre essentiel et le spectre singulier. Elle a également 
trouvé des applications étendues dans la physique, la théorie des équations aux dérivées 
partielles, l’analyse fonctionnelle et l'ingénierie. 

Aujourd’hui, la théorie spectrale est un domaine mathématique essentiel qui s’applique 
dans de nombreux domaines, allant de la physique à l’informatique quantique, en passant 
par l'ingénierie. Elle continue de se développer avec de nouvelles découvertes et applications. 
Elle joue un rôle crucial dans l’analyse de systèmes complexes et périodiques, en extrayant 


des informations sur les fréquences, les énergies et les modes de vibration. 
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1.2 Applications de la théorie spectrale 


La théorie spectrale trouve des applications diverses, telles que en mécanique quantique 
pour comprendre les opérateurs d’observable, en optique pour l’analyse de spectres lumi- 
neux, en ingénierie pour le traitement du signal et l’étude des ondes électromagnétiques, en 
informatique pour l’analyse de données et les algorithmes quantiques, en chimie pour la spec- 
troscopie, en économie et finance pour l’analyse des séries temporelles, et en biologie pour 
l’analyse de signaux biologiques. Le calcul du spectre d’un opérateur, en tant qu’application 
de la théorie spectrale, représente un puissant outil mathématique ayant des applications 
dans divers domaines. Il fournit des informations cruciales sur les propriétés spectrales des 
systèmes étudiés, jouant un rôle essentiel dans la compréhension, l’analyse, la modélisation 
et l’optimisation de systèmes complexes dans des domaines aussi variés que : 

Mathématiques : En analyse fonctionnelle, le calcul du spectre d’un opérateur est es- 
sentiel pour comprendre les propriétés des opérateurs linéaires, y compris leur compacité, 
leur inversibilité et la structure de leur spectre. 

Physique : En mécanique quantique, le calcul du spectre de l’opérateur hamiltonien est 
crucial pour déterminer les niveaux d’énergie possibles d’un système quantique, ce qui est 
essentiel pour comprendre le comportement des particules subatomiques. 

Ingénierie : Dans le domaine de l'ingénierie mécanique, le calcul du spectre d’un opé- 
rateur différentiel lié à un système mécanique permet d’analyser les modes de vibration, 
d'identifier les fréquences naturelles du système et de concevoir des systèmes anti-vibrations. 

Économie et finance : Dans la théorie des jeux, le calcul du spectre est utilisé pour 
analyser les stratégies et les équilibres dans des jeux à somme nulle, ce qui facilite la prise de 
décisions rationnelles dans des scénarios de négociation, d'enchères et de conflits. 

Géophysique : En sismologie, le calcul du spectre des opérateurs géophysiques permet 
de caractériser les propriétés des ondes sismiques, de détecter des structures géologiques 
souterraines, des tremblements de terre et des mouvements de plaques tectoniques. 

Informatique quantique : Le calcul du spectre des opérateurs est essentiel en informa- 
tique quantique pour la conception d’algorithmes quantiques, notamment pour la factorisa- 
tion d’entiers et la recherche quantique. 


Chimie : Dans le domaine de la spectroscopie, le calcul du spectre d’un opérateur quan- 
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tique permet de déterminer les transitions énergétiques et les niveaux d'énergie des molécules, 
ce qui est crucial pour l'identification de substances chimiques et l’étude des réactions chi- 
miques. 

Optique : En optique quantique, le calcul du spectre des opérateurs optiques est utilisé 
pour analyser la propagation de la lumière et la manipulation des états quantiques de la 
lumière. 

Biologie : En biophysique, le calcul du spectre d'opérateurs est employé pour analyser 
les propriétés des réseaux neuronaux et des systèmes biologiques complexes, ce qui peut 


contribuer à la compréhension du fonctionnement du cerveau et des cellules. 


1.3 Approximation du spectre et pollution spectrale 


Le calcul du spectre de la plupart des opérateurs non bornés, qui servent de modélisa- 
tion mathématique à des problèmes réels dans les domaines mentionnés précédemment, est 
souvent complexe et parfois presque impossible. C’est pourquoi les chercheurs ont recours à 
l’approximation du spectre pour obtenir des valeurs approximatives. 

L’approximation du spectre d’un opérateur non borné consiste à restreindre l’opérateur 
à un sous-espace de dimension finie, généralement en introduisant des conditions aux limites 
appropriées. Cette approche permet de travailler avec des opérateurs bornés, dont le spectre 
est plus facile à gérer mathématiquement. 

Cependant, l’utilisation de l’approximation du spectre d’un opérateur non borné peut 
entraîner un phénomène appelé pollution spectrale. Ce phénomène se produit lorsque l’ap- 
proximation inclut des valeurs propres qui ne font pas partie du spectre propre de l’opérateur 
non borné, désignées comme valeurs propres parasites ou artificielles. 

La pollution spectrale découle généralement des conditions aux limites ou des approxi- 
mations utilisées pour restreindre l’opérateur non borné. Des conditions aux limites inappro- 
priées ou des approximations grossières peuvent introduire des valeurs propres parasites qui 
ne correspondent pas aux propriétés physiques du système. 

Les techniques d’approximation du spectre d’un opérateur non borné varient en fonction 
du contexte. Certaines approches courantes incluent l’utilisation de séries de Fourier, de la 


théorie des perturbations, de la méthode des différences finies, de la méthode des éléments 
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finis et de la méthode de projection. 

Le phénomène de pollution spectrale constitue un problème récurrent dans l’approxima- 
tion du spectre d’un opérateur non borné. Il survient lorsque des approximations numériques, 
telles que des méthodes de discrétisation ou de projection, introduisent des valeurs propres 
fausses ou indésirables dans le spectre approximatif de l’opérateur, ce qui peut affecter la 


précision des calculs et la convergence des approximations numériques. 


1.4 Problèmes spectraux 


Les problèmes concrets dans les divers domaines mentionnés précédemment trouvent prin- 
cipalement leur solution grâce à l’étude des problèmes spectraux résultant de leur modéli- 
sation mathématique. Ces problèmes spectraux peuvent être regroupés en différentes caté- 
gories en fonction de leurs caractéristiques. Voici quelques-unes des catégories courantes de 
problèmes spectraux : 

Problème spectral standard : Il s’agit du problème de base où l'objectif est de trouver 
les valeurs propres et les vecteurs propres d’une matrice ou le spectre d’un opérateur linéaire 


donné. Cela consiste à résoudre l’équation suivante : 
Trouver (À,x), x £0, Ax = Àx 


où À peut être considéré comme une matrice, un opérateur borné ou un opérateur non borné. 

Ce problème a été traité dans le cadre d’une théorie spectrale classique, présentée en 
détail par Ahues et ses collaborateurs dans [I]. Dans ce document, les auteurs construisent 
une approximation de rang fini de À, transformant ainsi le problème initial Ax = Àx en un 
problème aux valeurs propres matricielle. Ils démontrent ensuite que la convergence de À, 
vers À en norme, en collectivement compact et en —convergence résout le problème de la 
pollution spectrale, établissant ainsi une approche numérique stable. 

Dans un autre contexte, les auteurs des références [IG], [IS], et [LA ont abordé le 
problème de la pollution spectrale qui survient dans la résolution d’un problème spectral 
standard associé à un opérateur non borné À. 

Dans ces travaux, les auteurs étudient un autre type de problème, à savoir 

le problème spectral généralisé, qui cherche à résoudre un problème spectral impli- 


quant simultanément deux opérateurs ou matrices. Les valeurs propres et les vecteurs propres 
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sont déterminés pour deux opérateurs liés, et le problème est formulé comme suit : 
Trouver (ic), z 20, Tr =)8x, 


où T'et S peuvent être considérés comme des matrices ou des opérateurs, qu’ils soient bornés 
ou non bornés. 
Ils démontrent que le spectre d’un opérateur non borné est égal au spectre généralisé de 


deux opérateurs bornés de la forme suivante : 
sp(A) = sp(T, 5), 


où, À est un opérateur non borné, T,S sont deux opérateurs bornés et les ensembles sp(A), 
sp(T, S) représentent respectivement le spectre de l’opérateur À et le spectre généralisé du 


couple (T°, S). Ces derniers sont donnés par 
sp(A) = {AE C/ À — ÀT n’est pas inversible }, 


et 


sp(T, $) = {1€ C/T — ÀS n’est pas inversible }. 


Les auteurs dans [I], [IG], [SI ont étudié le problème spectral généralisé Tx = ÀSx et ont 
démontré la convergence du spectre généralisé sp(T,,, $,) dans le cas de la convergence des 
suites 7, et $, vers T et S en norme [0], et en collectivement compacte [IGl, [8I. Tout cela 
est reconnu sous le nom de réalisation de la propriété U [I]. L'obtention de cette propriété 
montre que la pollution spectrale n’a pas de place dans ce type de problème. 

Le troisième type de problème est le Problème spectral quadratique, impliquant 
la résolution d'équations quadratiques en termes de valeurs spectrales. Ils sont définis par 


l'équation suivante : 
Trouver (À,x), x Z 0, XAx + ÀBzx + Cr =0, 


où À, B et C peuvent être des matrices, des opérateurs bornés ou des opérateurs non bornés. 
Ce type de problème occupe une place centrale dans la modélisation mathématique et 
est au coeur de nombreuses applications nécessitant la résolution de ce genre de problème, 


notamment dans le domaine de l’analyse dynamique des systèmes mécaniques et acoustiques, 
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la simulation de circuits électriques, la mécanique des fluides, ainsi que dans la modélisation 
de systèmes mécaniques microélectroniques. 

Ces problèmes ont attiré beaucoup d’attention dans divers domaines de la modélisation 
mathématique, tels que l’analyse dynamique des systèmes mécaniques en acoustique, la sta- 
bilité linéaire de l’écoulement en mécanique des fluides et la généralisation à la dimension 
infinie du problème du transistor (voir [I], [26]). 

Ces problèmes sont connus sous le nom de problèmes des valeurs propres quadratiques 
dans le cas où À, B et C sont des matrices, comme étudié dans [6[, [2], et [26]. 

Notons également que la résolution d’un problème des valeurs propres quadratiques joue 
un rôle très important dans la résolution d’une équation différentielle linéaire du second ordre 
telle que 

Mq"(t) + Cd'(t) + Ka(t) = (4), 


où, M,C' et K sont des matrices n x n et q(t) est un vecteur d’ordre n à déterminer. Cette 
équation illustre, dans le domaine des oscillations mécaniques, un système à ressort et masse 
unique (Figure 1) dans lequel un bloc rigide de masse M repose sur des rouleaux et ne peut se 
déplacer que par une simple translation. La résistance (statique) au déplacement est fournie 
par un ressort de raideur X, tandis que le mécanisme de perte d'énergie (dynamique) est 


représenté par un amortisseur C ; f(t) représente une force externe. 


FIGURE 1.1 — Système à masse unique-ressort. 


Dans le domaine de l'électricité, une équation différentielle linéaire du second ordre peut 
représenter l’écoulement du courant électrique dans un circuit RLC simple composé d’une 
bobine d’inductance L, d’une résistance À et d’un condensateur de capacité C, comme illustré 


dans la Figure 2; E(t) est la tension d’entrée. 


d?i(t) F RC ie dE(t) 


L 2 ol y 
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où i(t) est le courant à travers la résistance, et t est le temps écoulé. 


FIGURE 1.2 — Un circuit électrique simple. 


De même, dans un cas particulier, ces problèmes conduisent à l’étude des propriétés spec- 
trales du pencil quadratique de différents opérateurs non bornés, tels que le pencil quadratique 
de l’opérateur de Schrôdinger [BI et de l'opérateur de Sturm-Liouville [2], représenté par le 


problème différentiel suivant : 


Lay = —Y" + [q(x) + 2Ap(x)]y = Xy, (1.1) 
avec les conditions aux limites 


UD, À) = 1, d'(0, À) =D, 


ÿ{(r, À) + Hy(r, À) = 0, 


où, px), q(x) € L?[0, p]. 
Ce type de problème est connu dans la littérature sous le nom de problème de diffusion. 
De même, il décrit les interactions entre les particules lors d’une collision et présente un 


intérêt fondamental en physique [2]. 


1.5 Position du problème 


Dans notre travail, nous nous penchons sur le problème spectral quadratique associé à 


trois opérateurs linéaires bornés, formulé comme suit : 
Trouver (À,x), x £ 0, Ar + ÀBx + Cr = 0. 


Ici, À, B et C’ sont trois opérateurs linéaires bornés définis sur un espace de Banach B 


dans lui-même, avec À € C comme valeur spectrale. 
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Notre objectif principal est de démontrer l’absence de pollution spectrale dans l’approxi- 
mation du spectre quadratique généralisé de trois opérateurs bornés. Nous cherchons égale- 
ment à développer une nouvelle méthode numérique visant à minimiser et éviter la pollution 
spectrale lors de l’approximation de problèmes spectraux quadratiques, en particulier dans 
le cas des opérateurs non bornés. 

Considérons les suites d'opérateurs linéaires bornés 4,, B, et C, convergeant respective- 
ment vers À, B et C d’une certaine manière. Deux questions essentielles se posent naturelle- 
ment : 

— Si À E SP(Ay, Ba, Cn) et An — À, est-ce que À € sp(A, B,C)? 

— Sie sp(4, B,C), existe-t-il À, € sp(A,, B,, C,) pour chaque grand n tel que À, — À? 

La réponse négative à la première question définit la notion de la pollution spectrale. En 
revanche, la réponse négative à la deuxième question montre que le spectre n’atteint pas sa 
totalité. 

Notre objectif à travers cette thèse est de rechercher les modes de convergence de 4,, B, 
et C, vers À,B et C’, qui permettent d'obtenir une réponse positive à ces questions. Cela 
nous permettrait de conclure que les deux propriétés suivantes sont réalisées, résolvant ainsi 
le phénomène de la pollution spectrale et montrant que le spectre sera atteint dans sa totalité 


dans notre cas : 


La propriété U est vérifiée si À,, B, et C, convergent respectivement vers À, B et C dans 


un sens donné, À, € sp(A», B,, Ch) et An — À. Alors À € sp(A, B,C\). 


La propriété L est vérifiée si 4,, B, et C, convergent respectivement vers À, B et C dans 
un sens donné, et que À € sp(À, B,C). Il existe un certain À, € sp(A,, B;,C») pour 


chaque n suffisamment grand, de telle sorte que À, — À. 


1.6 Plan de la thèse 


Notre thèse se divise en trois chapitres essentiels. Le premier chapitre sert de fondement et 
introduit des concepts préliminaires importants pour la compréhension des chapitres suivants. 
Ces concepts comprennent un rappel du spectre généralisé et des résultats obtenus lors de 
l’approximation du spectre généralisé. De plus, nous présentons quelques notions et résultats 


relatifs au spectre quadratique généralisé, établissant ainsi le cadre de notre étude dans les 
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chapitres ultérieurs. 

Le deuxième chapitre est consacré à la présentation de la méthode d’approximation du 
spectre quadratique généralisé. Nous examinons la démonstration des propriétés U et L, où 
les suites d'opérateurs À,, B, et C, convergent en norme et collectivement compacte vers les 
opérateurs bornés À, B et C, respectivement. Un élément crucial de notre travail consiste à 
développer une extension de la notion de convergence classique v pour démontrer la propriété 
U associée aux opérateurs À, B, et I + D. 

Dans le dernier chapitre, nous appliquons la méthode du spectre quadratique généralisé au 
problème de l’opérateur de Schrôdinger sous forme de pencil quadratique, afin de résoudre le 
problème de pollution spectrale qui se pose lors de l’approximation de ce type d'opérateurs 
non bornés. De plus, nous réalisons des tests numériques pour valider l'efficacité de notre 


méthode. 
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CHAPITRE 2 


Spectre généralisé et spectre 
quadratique généralisé 


SomiA. KAMOUCHE UNIVERSITÉ 08 Maï 1945 GUELMA 


2.1 Introduction 


Ce chapitre est consacré à l’introduction des fondements théoriques utilisés dans le trai- 
tement du spectre quadratique généralisé, associé aux trois opérateurs linéaires bornés, ce 
qui va nous permettre de developper une nouvelle méthode d’approximation du spectre qua- 
dratique généralisé. Aussi, ils servent à assurer la convergence de cette nouvelle méthode 


moyennant des hypothèses appropriées. 


Ce chapitre est organisé selon le schéma suivant : 


Un rappel sur les résultats obtenus dans [9], concernant le spectre généralisé associé aux 
opérateurs bornés ; puis nous présentons une étude réalisée sur le spectre généralisé avec un 


nouveau mode de convergence appelé v—convergence généralisée. 


Ensuite, nous introduisons les nouvelles notions du spectre associé aux trois opérateurs 


linéaires bornés, où ces notions sont déjà initialement définies pour trois matrices [26]. 


Enfin, nous donnons les propriétés du spectre quadratique généralisé et de la résolvante 
quadratique généralisé, où ces propriétés sont la généralisations des concepts théoriques des 
cas classique et généralisé [I], [O, HI, [Gl, El, HA puis pour terminer avec une étude 
de la différentiabilité de la résolvante quadratique généralisé au sens de Fréchet qui va nous 


permettre de montrer des propriétés plus régulières de l’opérateur résolvant par rapport à 


celles obtenues dans [E3l, [A], [F1]. 


2.2 Rappel de quelques notions et résultats sur le spectre 
généralisé 


Soit (B,|| -|}#) un espace de Banach. On considère l’espace noté BL(B) qui est l’ensemble 


des opérateurs linéaires bornés de B dans B, i.e. 


BL(B) = {T:B-—8B, T linéaire borné}. (2.1) 
On munit cet espace de la norme notée || - || définie, pour T € BL(B), par 
IT = sup{ÎTuls : u € B, lus = 13. (2.2) 
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Définition 2.2.1 


Soient T et S deux opérateurs dans BL(B). 


1. On appelle problème du spectre généralisé, le problème donné par 


Trouver À € C etu € B—{0} tel que Tu = ÀSu. (2.3) 


2. On appelle spectre généralisé des opérateurs T et S, l’ensemble noté sp(T, S) défini par 
sp(T,S)={1EC/T —AS n’est pas inversible }. (2.4) 


3. On appelle ensemble résolvant généralisé des opérateurs T et S l’ensemble notére(T,S), 
défini par 
re(T, S) =C\sp(T,S). (2.5) 


4. On appelle spectre ponctuel généralisé des opérateurs T et S, l’ensemble noté sp,(T, S) 


défini par 


sp,(T, S) = {A € C/3u € B\{0}, Tu = ÀSu}. (2.6) 


5. On appelle opérateur résolvant généralisé de (T,S) en z noté R(z,T,S), la fonction 


défini par 


R(;,T,S):re(T, S)CC — BL(B) 


2 + RGTS=(T-3238) 7 


6. On appelle élément spectral généralisé isolé, tout élément À dans sp(T, S), tel qu'il existe 


une courbe dans re(T, S) notée T et appelée contour de Cauchy, et 
Xe int(T), (2.7) 


avec la condition 


int(T) Nsp(T, S)\ {À} = ÿ, (2.8) 


où int(l) désigne l’intérieur du contour de Cauchy T. 
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Les deux théorèmes suivants servent à introduire les outils fondamentaux, pour l’étude 
du spectre généralisé et la mise en relief de quelques propriétés essentielles, à savoir que 


l’ensemble du spectre généralisé est fermé et que la résolvante généralisé est analytique. 


Théorème 2.2.1 
Soient T et S deux opérateurs linéaires bornés, T, S € BL(B). On a alors les résultats sui- 
vants : 
1. Soit z appartenant à l’ensemble résolvant généralisé re(T, S) et soit z0 appartenant à 
C. 
Si|z — zol < IIR(2,T, SSI, (2.9) 
alors, 20 € re(T, S). 


2. Le spectre généralisé sp(T,, S) est un ensemble fermé de C. 


3. La fonction résolvante R(-.,T,S) : re(T, S) — BL(B) est analytique et sa dérivée par 


rapport à 2 vérifie la relation : 


d 
TG T,5)=R(2T,SYSR(AT,S). (2.10) 


4. Soient T, S,T,S dans BL(B), et z € re(T, S). 


Si les deux inégalités (2.11) et (2.12) ci-dessous sont vérifiées 
ICT) - 28 — 5) R(2,T, 5) < 1, (2.11) 
où 


L(T-n 248 -8))R(2T,8)] | <1, (2.12) 


alors 2 est dans l’ensemble résolvant de T et S (2€ re(T, S)). 


Preuve. Voir 9]. 


Théorème 2.2.2 


Considérons deux opérateurs bornés T et S, T,S € BL(B). Soit B(0,k) la boule de centre 0 
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et de rayon k > 0. 


On a alors les assertions suivantes : 


1. L'ensemble du spectre généralisé sp(T, S) est un sous- ensemble de la boule B(0,k), si 


et seulement si, 0 n’est pas un élément du spectre de S, 1.e. 
sp(T, S) € B(0,k) & 0 é sp(sS). (2,13) 
2. Si l'opérateur S est compact, alors le spectre généralisé sp(T,, S) est la réunion du spectre 
ponctuel généralisé sp,(T, S) et l'infini, 1.e. 
sp(T”, S) = sp, (T°, $) U {oo}. (2.14) 


Preuve. Voir A9]. 


Le théorème suivant carctérise l’existence de l’opérateur projection. 


Théorème 2.2.3 
Soit À une valeur propre généralisée de type fini et isolée, dans le spectre généralisé sp(T, S). 


Si T' désigne le contour de Cauchy qui sépare À de sp(T, S), alors l'opérateur 


—] 
= — T = 28) Sd 2.1 
est une projection, ï.e. P? = P. 


Preuve. Voir [1]. 


2.3 Approximation du spectre généralisé 


Dans cette section, nous présentons un schéma pour l’approximation du spectre généralisé, 
Ensuite, nous énonçons les principaux résultats, qui assurent la convergence de la nouvelle 


méthode du spectre généralisé, sous des hypothèses appropriées. 


Soient (Tn)nen et (Sn)nen deux suites d'opérateurs linéaires bornés, qui convergent dans 
un sens précis vers les opérateurs bornés T' et S respectivement. Alors, pour déterminer la 
relation entre sp(T, S) et sp(T,, S,) (pour n assez grand), nous considérons le schéma d’ap- 


proximation spectrale élaboré à travers les propriétés suivantes : 
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Propriété U : Si À, € sp(Th, Sn) et À, — À dans C, alors À € sp(T,S). 


Propriété L : Si À € sp(T,S), alors il existe une suite (X,),en telle que À, — À dans 
C et }n € SP(Tn, Sn). 


Remarque 2.3.1 


1. Notons que la Propriété U et la Propriété L sont les extensions naturelles du cas stan- 


dard S = I, T € BL(B) et I est l'opérateur identité dans BL(B) (voir O, HW). 


2. Pour clarifier la notion du "sens précis" de la convergence, nous citons quelques modes 
de convergence usuelles : 


— La convergence ponctuelle, notée Ty + T, telle que 


Pour toutz € B, | Tax — Tx | 0. 


/ nm 
— La convergence en norme, notée T, — T, telle que 


ÎTn —T | 0. 


° a CC Ze 
— La convergence en collectivement compact, notée T, = T, vérifiant 


TT 
et pour un certain nombre entier positif no 


U {Tr-Tr ire, (els <1} 


n>no 


est relativement compact dans B. 
— La v — convergence, notée T, + T, définie par 


sup [7h] < +oo, et on a les convergences en norme 
nEN 
(Tr — TT] — 0 et (Ts — TT] — 0. 
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Dans toute la suite, on montre que le schéma de la convergence du spectre généralisé est 
obtenu à l’aide des propriétés U et L, en faisant appel à la convergence en norme et en 


collectivement compact. 


Théorème 2.3.1 
Soient {Th }nen Et {Sn}nen deux suites d'opérateurs bornés. Si les opérateurs T,, et S, vérifient 


l’une des hypothèses suivantes : 
(A1) T, ST, Sn “+ S. (convergence en norme) 
(A2) T, ST, S, + S. (convergence en collectivement compact) 


et si pour tout n E N, Ày E SP(Tn, Sn) et À, converge vers À dans C, alors À € sp(T, S). 


Preuve. Voir A9]. 


Théorème 2.3.2 

Si les suites d'opérateurs {Ti }nen et {Sn}nen vérifient l’une des hypothèses suivantes : 
(A1) M ST, Sn +S$. 
(A2) T, ST, Sn + S. 

et si À est une valeur propre généralisée isolée de type fini dans sp(T, S), alors il existe une 


suite ]n € SP(Tn, Sn) telle que À, converge vers À dans C. 


Preuve. Voir H9]. 


Dans le but d'améliorer et d’affaiblir les hypothèses du théorème et du théorème 
E.3.2 nous allons introduire un nouveau type de convergence plus faible que les autres modes 
de convergences utilisés dans HI], [A. 

Nous appelons ce type de convergence la y — convergence généralisée associée aux opérateurs 
T, et 5, . Cette v — convergence généralisée est définie de la manière suivante : 

On dit que (7,,$5,) converge vers (T,S) en v — convergence généralisée et on note 

(Tu, Sn) + (T, S) si les conditions suivantes sont vérifiées : 


— |T,ll et |S,1] sont bornés pour tout n, on a alors 


VneEN, [Tall < æet |[S,| < oo. (2.16) 
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— Les quantités suivantes doient être tendre vers 0 lorsque n — oo 
nr — T) Tai = En — T) Sal = (Sn — 5) Snll = (Sn — S)Tnl +0 (2:17) 
In = TT = (Tr — TT) SÙ = (Sn — S) 51 = 1 (Sn — SIT > 0: (2.18) 
Dans toute la suite, nous allons traiter le spectre généralisé sp(T, S), avec T=1+@G 


où J est l'identité, $,G sont deux opérateurs bornés, $,G € BL(B). De même, le problème 


spectrale généralisé (2.3) s'écrit 


Trouver À € C et u € B\{0} tel que (1 + Gju = ÀSu. (2.19) 


Le problème (2.19) sera traité en adoptant la y — convergence généralisée. 


Remarque 2.3.2 
Les résultats de cette nouvelle approche du problème du spectre généralisée à l’aide de la 


y — convergence généralisé a fait l’objet d’une publication, (voir[11)). 


Les résultats essentiels du traitement de notre problème spectrale généralisé (2.19) sont 
illustrés à travers les deux théorèmes ci dessous, où la propriété U a été prouvée sous la 


v — convergence généralisée. 


Proposition 2.3.1 
Soient G et S dans BL(B). 


(a) Pour z1 et z2 appartenant à l’ensemble résolvant re(1 + G,S), on a alors 
R(4) — R(2) = (4 — 2)R(4)S R(2). (2.20) 
(b) Pour z € re(l + G,S), on a alors 
R(2)=1-GR(:) +2SR(). (2.21) 


(c) S'il existe z0 € re(1 + G,S), tel que l'opérateur G commute avec l'opérateur resolvant 


généralisé au point z0 R(20),t.e. 
alors 
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(c.1) Pour tout z et 22 € re(l + G,S), l'opérateur R(z1) commute avec l'opérateur 


R(2), R(2),1.e.Rp(z1) = R(a) R(2). 


(c.2) Pour tout z E re(1+G,S), les opérateurs G,S commutent avec l'opérateur R(2). 


Preuve. 


Soient S et G deux opérateurs dans BL(B). 
(a) Soient 21,22 deux élément de l’ensemble résolvant quadratique généralisé re(1 + G,S). 
On écrit 
R)=R&)=U1+G- AS) T-U1+G-328)" (2.23) 
On met R(z) comme facteur commun, pour obtenir 
R(a) —-R(z2) = R(4)1-(1+G-24S)R(2)]: 


De plus, on met R(22) comme facteur commun, et on obtient 


R(ä) Rx) = R(4)[1+G-28)-(1+G-—2S)R(2) 


(z1 — 22) R(41)S R(2). 


Ceci montre le cas (a) du théorème. 


(b) Pour z € re(J + G,$), on peut écrire l’opérateur identité par la formule suivante : 
[1 = (1+G—-2S)R(:) =R(:) +GR(:) — 2SR(2), 
ce qui permet d'écrire l’opérateur R(z) comme suit 
R(2)=1-GR(:) +2SR(:). (2.24) 


Cela prouve (b) du théorème. 
(c) Soit & Ere(1+G,S). 


(c.1) Soient 2, et z, deux éléments de re(J + G, 5). Comme R(z1) commute avec S et 


G, alors, on peut écrire 


R(z1) [7 + G— 225] — [7 + G — 225] R(21) 
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En multipliant par R(22) à droite et à gauche, on obtient 


R(22) R(41) = R() R(2). 


Par conséquent, (c.1) est vérifiée. 


(c.2) En utilisant le résultat obtenu dans (a), pour tout z € re(J + G,S), et pour 


zo € re(1 + G,S), on montre que les opérateurs $S et G commutent avec R(2). 


Le théorème suivant montre que la proprieté U est vérifiée pour la v — convergence géné- 


ralisée. 


Théorème 2.3.3 
Soient Gn et S, deux suites d'opérateurs bornés, et À, un élément du spectre généralisé 
Sp(l + Gn, Sn) , telle que le couple (G»,$,) converge selon la v — convergence généralisée, et 


An — À dans C. Alors À est un élément de sp(1 + G,S). 


Preuve. 
La démonstration de ce théorème est basée sur le raisonnement par l’absurde. On suppose 
que À € re(1 +G,S). Comme l’ensemble résolvant généralisé re(J + G, S) est un ouvert dans 


C, alors il existe un scalaire Ô > 0 tel que 
Q={z2eC:lz2-AÀ|<ô} Cre(l+G,S). 


De plus, la résolvante généralisée R(-) étant analytique sur l’ensemble re(J + @,S), on peut 


définir le nombre réel 5 par 
a = sup{||R(2)|| : z € Q} < oo. 


On cherche à montrer que l’ensemble Q est inclus dans l’ensemble re(]+@G,,$,). Pour cela, il 
suffit de prouver que pour tout z € A on a z € re(1 + Gh, Sn), i.e. on montre que l’opérateur 


Ra(z) = (1 + Gn — 2S,) ! existe et est borné, pour tout z € (1. 
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D'une part, on va développer le terme [[(G — G,) — 2(S — S,)]R(z:)|, en utilisant l’in- 
égalité présentée dans la proposition b), qui nous permet d'écrire R(z) sous la forme 
suivante : 


R(2) = 1-GR(z:) +2SR(:). 


Cela donne 


[CG — Ga) — 28 — SR) 
= [(G-G)-2(5- SIT -GRQ() +2SR()II(G — Gr) — (5 — S)IR() 


LS Ra(z) — GRa(2)][(G — Gr) — 2(S — S,)IR(z) 
et comme |G,|| et ||S,|[sont bornés, ils existent t, s > 0 tels que |G—G,|| <tet |S—S,] < s 
pour n € N. D’après la convergence des opérateurs G,; et S, en v—convergence généralisée, 


on à 


I (Gn — G)Gnll = (Sn — 5) Gall = (Gr — G) Sn = 1 (Sn — 9) Sn — 0, 


et 


I (Gn — G)G| = (Sn —S) GT = (Ga — G)S1 = (Sn — S)S1 > 0. 


Pour montrer que RA(2) = (1 + Gn — 2$,) ! existe et est borné, on doit établir que 


AG) - 25 - SAR] < 1. 
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Par conséquent, pour tout z € (, on suppose qu'ils existent des entiers naturels n;, 1 < i < 8 


et & = max(t,s) vérifiant les hypothèses suivantes : 


1 
Ie EG TC EST NUE 
WG G)ST < Croate 
1 
EST EE Greta) 2" 
1 
MES E Gorspror sep "274 
I(G, - GG] Sy 2% 
(E-G)S < rEnp 27e 
(SSSR "277 
(SSD < ppp "Er 


Ceci nous permet de vérifier que 
JG — G) - 25 - SRG] < 1, 
et on à, pour tout n > max{n;,1 <i<8}, 
Ie -G)-28- SAR < 5 <2 


Alors 
Ra(z) = +Gn— 25) 


existe et est borné; donc z € re(1 + G,,$,). Par conséquent, Q C re(1 + G,,S,) pour tout 
n assez grand. 

Comme la suite des valeurs spectrales approchées À, converge vers la valeur exacte À, on à 
Ân EN Cre(l +G»,$Sh) pour tout n. On à aussi À, € sp(l + Gh,$,) pour tout n, et donc 
À doit être dans sp(J + G,S), ce qui est en contradiction avec l’hypothèse de départ que À 
n’est pas dans sp(1 + G,S). 
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2.4 Quelques notions et résultats sur le spectre qua- 
dratique généralisé 


Dans cette section, on va introduire des nouvelles notions sur le spectre quadratique gé- 
néralisé et des outils qui seront utilisés dans la construction de la nouvelle méthode d’ap- 
proximation de ce dernier. Cela permet l’extention de l’étude du spectre généralisé au spectre 
quadratique généralisé, ainsi qu’une généralisation des concepts connus dans les problèmes 


de valeurs propres quadratiques. 


Définition 2.4.1 


On considère trois opérateurs À, B et C' dans BL(B). 


1. On appelle problème spectral quadratique généralisé, noté Q()), l'équation 
Trouver À E C etu € B\{0}, Q(A)u = 0, (2,25) 


où, 
Q()u = X'Au + Bu + Cu. (2.26) 


2. On appelle ensemble résolvant quadratique généralisé des opérateurs À, B et C’, l’en- 


semble noté re(A, B,C') défini par 
re(À, B,C) = {2 EC: QX) = (A+2B +0C) est inversible et borné} . (2.27) 


3. On appelle opérateur résolvant quadratique généralisé du triplet (A, B,C) en z, noté 


Ra(z, À, B,C), la fonction définie par 
Ra(-) = Ra(:,4,B,C):re(A,B,C)CC — BL(B) 
z + Ro(z) = [QU] = (224 + 2B + CC). 
4. On appelle spectre quadratique généralisé associé aux opérateurs À, B et C notésp(A, B,C), 
le complémentaire de l’ensemble résolvant quadratique généralisé dans € donné par 
sp(À, B,C) = C\re(A, B,C\). (2.28) 


5. On appelle spectre ponctuel quadratique généralisé des opérateurs À, B et C', l’ensemble 


noté sp, (A, B,C')et défini par 


sp, (4, B,C)= {\E€ C, a € B\{0}: Q(Vu = VAu+ABu+Cu=0}. (2.29) 
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6. On appelle une valeur spectrale quadratique généralisée isolée toute valeur À dans sp(A, B, C) 
telle qu'il existe une courbe dans re(A, B,C) notée T, appelée contour de Cauchy, qui 
sépare À de sp(A, B,C')\ {À}, et 

Xe int(T) (2.30) 
avec 


int(T) Nsp(A,B,C) \ {A} = 6, (2.31) 
où int(l) désigne l’intèrieur du contour de Cauchy T. 


7. On appelle rayon de convergence le nombre noté K(À) défini par 


IR ICAIAIIBD) +41 ReCOIIAI 
HA) 2 RaQ IA F8 
|RaQICIAIAI + 181) 


2] Re() INA 


pour tout À € re(A,B,C. 


2.5 Propriétés de la résolvante quadratique généralisée 


Nous allons procéder à l’extension et la validation des propriétés du spectre classique et 


du spectre généralisé à notre spectre quadratique généralisé. 


Le théorème suivant montre que, même s’il y aura une perturbation d’un élément de 


l’ensemble résolvant quadratique généralisé, on a toujours une certaine stabilité de ce dernier. 


Théorème 2.5.1 
Soient À, B et C' trois opérateurs dans BL(B), soit z0 € C. Si z est dans l’ensemble résolvant 


quadratique généralisé re(A, B,C) tel que 
2 — 20 < k(2) 


où K(z) est donné par la formule E£.5, alors, 2 € re(A, B,C). 


Pour démontrer ce dernier théorème, on a besoin d’énoncer le théorème suivant qui sera 


aussi utilisé dans le deuxième chapitre. 
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Théorème 2.5.2 


Soit T un opérateur borné dans BL(B). 
Si la norme ||T|| < 1 alors (T — I)? existe et est borné. (2.33) 


De même, on a la formule suivante (série de Neumann) 


(T-1'=-S IT" (2.34) 


n>0 


et la norme ||[(T — 1)-!|| est donnée par 


1 


(TD Tr 


(2.35) 


Preuve. 
Pour montrer que z € re(A, B,C), il suffit de montrer que l’opérateur résolvant quadratique 
généralisé Ra(20) existe et est borné. 


Soit l'opérateur Q(z) donné par 
Q(20) = HA + 2%B +C. (2.36) 


Pour montrer l'existence de l'opérateur RQ(20), il suffit de montrer que l'inverse [Q(20)]7! de 


l'opérateur Q(2) existe. D'abord, on écrit 


Q(20) = Q(2) — Q(z) + Q(:) (2.37) 
(2.38) 

et d’après (2.37), on a 
Q() = ($-2)A+(x-2B+Q(2 (2.39) 


On met Q(z) en facteur dans (2.39) pour obtenir 
Q(z0) = Q(2) [1 — [Q(2) "(27 — 2) A + (2 — 20)B)| (2.40) 
et comme z € re(AÀ, B,C'), ceci donne 


Q(&0) = Q(2) [I — Ra(z) ((27 — 2)A + (2 20)B)] (2.41) 
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D'après (2.41), l'inverse de Q(20) existe si l'inverse de l’opérateur 
LT — Ro(z) (2? - 22)A + (2 - 2)B)| (2.42) 


existe. 
Prouvons que l'inverse de ce dernier existe. En utilisant la condition |z2 — z9| < «(z), on 
obtient 

I Ra(z) (27 — 22) + (2 — 20)B) || < 1. (2.43) 


Selon (2.33) du théorème [2.5.2| l'opérateur 
M 25 sn 
[LT — Ro(z) (2 - 2) + (2 - 2)B)| (2.44) 
existe et est borné. Par conséquent sa norme est majorée par 


[1 RoU) (2 - A+ 0 X0)8)] 7 | (2.45) 


< (1—IRaG) (A -DA+ (A X)B) 1) 


Comme z € re(A, B,C), l'opérateur résolvant quadratique généralisé RQ(z) existe et est 


borné; ce qui implique que 
[Q(z0)] ? = Ro(zo) = [1 - Ra(z) (2 - 2)A+(2- 2)B)]  Ra(). (2.46) 

De (2.45), la norme de Ro(z) vérifie 
IRa(zo)ll € (1— II Ra(z) (22-22) + (2 — 20)B) 1) IRa(z)ll (2.47) 


Donc, RQ(%) existe et est borné, ce qui montre que à € re(A, B,C). 


Corollaire 2.5.1 


Le spectre quadratique généralisé sp(A, B,C) est un ensemble fermé dans C. 


Preuve. 
En utilisant le fait que re(A, B,C') est le complémentaire de sp(A, B,C') dans C, et d’après 


le théorème précédent, nous obtenons que re(A, B, C') est ouvert dans C. 
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Théorème 2.5.3 

L'opérateur résolvant quadratique généralisé RQ(-) est une fonction analytique sur re(A, B,C). 
De plus, sa dérivée par rapport à z est donnée par l'expression suivante : 

pour tout z e re(A, B,C), 


dRa 
dz 


(z) = —Ro(z)(22A + B)Ra(z). (2.48) 


Preuve. 
Soient z, z0 deux éléments de l’ensemble resolvant quadratique généralisé re(A, B, C') tels que 


[z — 20l < k(z0). On écrit 
à di $ =i 
Ra(z) —Ro(z) = (24+:B+C) -(#4+2%B+C) . (2.49) 
On prend RQ(z) comme un facteur commun à droite et RQ(20) à gauche, pour obtenir 
Ra(z) — Ra(zo) = Ra(z) ((2ê — 2?) A + (20 — 2)B) Ra(zo). (2.50) 


Pour montrer que l’opérateur RQ(z) est analytique et sa dérivée est donnée par (2.48), il 


suffit d'établir que 


lim Rat — _ E Ro(z)(2:A + B)Ro(z)| = 0. (2.51) 
On commence à développer le terme 
| Ra() = Ra) | R (211224 + BYRo(e) 
Z — 20 
En utilisant la relation suivante : 
Ra(z) — Ra(zo) = Ra) ((4 — ?)A + (20 — 2)B) Ro(zo), (2.52) 
on obtient 
| Rat — _ LRa(z)(2:A + B)Ro(z) 
: | LONGER ONG) Era), CS 
— 20 
Aprés simplification par z — 20, on obtient 
| Rat) L ae L Ra(z)(2zA + B)Ra(z) 
= Î-Ra()((o + z)4 + B) Ra() + Ra(z)(2:A + B)Ra(z)| (2.54) 
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Ensuite, on ajoute et on retranche des termes, puis on applique les propriétés de la norme. 
Cela permet d’obtenir l’inégalité suivante : 


Ra(z) — Ra(o) 


2 — 20 


zo + 2] ]Ra(z)AÏ IRa(z) — Ra(zo)| + IRa(z)BITIRa(z) — Ra(zo)| 


+ 12-201] Ra(z)A Ra(2)|| 


+ Ro(z)(2:A + B)Ra(z) (2.55) 


PA 


Pour établir que le terme 


| Ra(z) — Ra(z0) 


2 — 20 


LE Ra(z)(2zA + B) Ra(z) 


tend vers 0 quand z tend vers 29, il suffit de montrer que 
IRa(z) — Ro(zo)| 
tend vers 0 aussi quand 2 tend vers 20; pour cela, en utilisant la formule (2.44) 
2.5 a 
Ra(zo) = [7 —Ra(z) ((2?—28)4+(2- 2)B)] Ro(z) 


obtenue dans la preuve du théorème |2.5.1| et à l’aide du dévelopement en série de Neumann 


(2.34) du théorème donné par 


[1 — Ro(z) (2 - 2)A + (2- 2)B)| = Ro)" (22-22) A+ (2 20)B) . (2.56) 


On a alors, h 
IRa(z) — Ro(z)|| = (prit "(22 - 2)A+ (2 2)B) Ro) 
_ - Z Rat) 1 - a+ 208)" Ref) 
= - Z Rat) "(2 — 20)" ((z+ 20) A + B)" Ra(z) 
< XL IRaG)( — 20) (2 + 20) + BR GI. 
Comme 


IRa(z)(z — 20) ((z + 20) À + B)] <1, 
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cette série géométrique converge, ce qui nous donne 


_ IRa(2)II IC — 20)! (Ce + 20)4 + BI IRo(2)I 
IRoG) RG ET I + 2094 + BI IR GI 


Cette dernière inégalité montre que 
Jim [Ra(z) — Ra(z)|| = 0. (2.57) 


Cela implique que 


Ra(z) — Ra(2o) 


Z — 20 


lim 
Z—20 


Ra(z)(2zA + B) Ra(z) 


= 0. (2.58) 


Par conséquent la dérivée de l’opérateur RQ(z) par rapport à z est 


Te (2) = —Ro(z)(2zA + B)Ro(2). 


Le théorème suivant est un résultat essentiel qui nous permet de montrer la propriété U 


sous notre nouveau type de convergence (la 7—convergence généralisée). 


L'idée fondamentale est d'écrire l'opérateur C = 1 + D dans (2.26), où D est un opérateur 
borné dans BL(B). 


Théorème 2.5.4 


Soient À, B et D trois opérateurs dans BL(B). Les propriétés suivantes sont vérifiées 


(a) Pour z1 et z2 dans l’ensemble résolvant quadratique généralisé re(A, B,1 + D), on a 
Ra(z) — Ro(z2) = Ra(a) [(23—)A+(2-2)B]Ra(z). (259) 
(b) Pour z € re(A, B, I + D), l'opérateur RQ(2) vérifie 
Ra(z) = 1 —- 2 ARa(2) — zB Ra(z) — DRa(z). (2.60) 


Preuve. 


Soient À, B et D trois opérateurs dans BL(B). 
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(a) Soient z1, 22 deux éléments de l’ensemble résolvant quadratique généralisé re(A, B, 1+D). 


On écrit 
Ra(z) — Ro(z) = (A+4B+1+D) -($A+2B+1I+ D). (2.61) 
On prend RQ) comme facteur commun, pour obtenir 
Ra(a)—Ro(z) = Ro(a)[1-(#4+248 +14 D)Ro(x)]. 
De même, en mettant RQ(2) en facteur, il s'ensuit 


Ra(z1) — Ra(z2) 


Ro(a) [(24+2B+1+D)-(#A+4B+1+ D) Ro(z) 


= Ro(z) |(23 - 2)A + (2 - 4)B] Ra(z). 


Ceci montre le cas (a). 


(b) Pour z € re(A, B,1 + D), on peut écrire l’opérateur identité par la formule suivante : 
I = (A+2B+1I+D)Ro(z) = ARo(z) + zB Ro(z) + D Ro(z) + Ro(z). 
Alors, l'opérateur RQ(z) s'écrit 
Ra(z) = 1-2 ARa(z) — 2BRa(z) — DRa(2), (2.62) 


ce qui prouve (b). 


Le théorème suivant représente un résultat essentiel qui nous permet de montrer la pro- 


priété L du schéma d’approximation du spectre quadratique généralisé. 


Théorème 2.5.5 


Soient À, B et C' trois opérateurs dans BL(B). 


(a) Pour z € re(A, B,C') \ {0}, on a alors 
(B + 2A) Ra(z)C = CRa(z)(B + 2A). (2.63) 


(b) Si l'opérateur C' commute avec RQ(z1) et RQ(22), pour certains 1,22 € re(A, B,C)\{0}, 


alors C commute avec À, B et RQ(z) pour tout z € re(A, B,C). 
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Preuve. 


(a) Pour tout z € re(A, B,C') \ {0} on va montrer la relation 
(B +2A)Ra(z)C = CRa(z)(B +24), 


en développant ses deux membres . 


Pour le membre de gauche, on prend 1/2 comme facteur commun pour avoir 
1 
(B+zA)Ro(z)C = -(2B+zA)Ro(z)C, (2.64) 


2 


Ensuite, on ajoute et on retranche le même terme pour obtenir 


ÉaRo= (4 + 2B + C) Ro(2)C — LC Ra(2)C (2.65) 


1 1 


De manière similaire, on développe le membre de droite. On obtient 


CRal(2(B+z4) = 20 Ro(z)(zB + 24) (2.66) 


1 1 
= -CRa(z)(7A +zB+C)— = CRa(z)C 
1 1 
De (2.65) and (2.66), on a 
(B + zA) Ra(z)C = CRa(z)(B + zA), 


ce qui montre le cas (a) 


(b) Pour tout z; € re(A, B,C) \ {0},i= 1,2, on a C'Ra(z) = Ra(x)C, ie. 
C(ÉA+2B+C) "= (ZA +2B +C)C. 
Multiplions à gauche et à droite par (z?A + 2,B + C), pour obtenir 
(22CA + zCB + C?) = (AC + z,BC + C?), 
ce qui permet d'écrire 
4 [4(CA — AC) + (CB — BC)] = 0. 
Par conséquent, AC = C'A et C'B = BC: 
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Remarque 2.5.1 


1. Le résultat obtenu en (a) est une généralisation de certains résultats présentés par 


Gohberg [A]. 


2. L'hypothèse que l'opérateur C commute avec les opérateurs résolvants quadratiques gé- 
néralisés RQ(4) et Ra(z2) pour certains 21,22 dans l’ensemble re(A, B,C') \ {0}, peut 
être utilisée pour l’opérateur À au lieu de l’opérateur €, ce qui donne le même résultat. 
En fait, on peut réécrire le problème spectral quadratique 2? A + 2B + C' sous la forme 
Q(2) = A+23B + 220 = 0 avec À = 1/2, 2 £ 0. 

3. Les deux derniers théorèmes à savoir théorème|2.5.1 et théorème |2.5.ë] représentent les 


outils principaux utilisés pour développer la deuxième section du chapitre 2. 


2.6 Dérivée de Fréchet de la résolvante quadratique gé- 
néralisée 


Rappelons que la méthode de Newton-Kantorovitch est une méthode bien connue dans 
le calcul approché des valeurs propres matricielles, puisque qu’ elle permet également le 
développement des méthodes du raffinage itératives des valeurs propres approximatives. Pour 
cette raison, nous pensons à utiliser cette méthode pour faire le raffinage des valeurs propres 
quadratiques généralisées. Pour cela, nous aurons recours à l’étude de la dérivée de Fréchet 
de l’opérateur résolvant quadratique généralisé qui est considéré comme une fonction non 
linéaire sur BL(B). Cette étude nous permet de montrer des propriétés plus régulières de 


l'opérateur résolvant par rapport à celles obtenues dans [T3]. 


Pour introduire la notion de la dérivée de Fréchet, on considère deux espaces de Banach, 
notés B1 et BB. 
On note par BL(B:1, B2) l’espace des applications linéaires bornées, définies de B; dans B2. 


Soit F une application définie sur Q, un ouvert de B, telle que 


F:0C B:1 — Bo. 
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Définition 2.6.1 
On appelle dérivée au sens de Fréchet de F l'opérateur linéaire continue L € BL(B:, B:) 
défini de la manière suivante : 


Pour x € Q et pour tout h € Bi. 
L(x)h = F(x+h) —F(x) —o(|h|). 


où o(||h|]) tend vers O quand ||h|| tend vers 0. 


Le théorème suivant nous permet le calcul de la dérivée de Fréchet de l’opérateur résolvant 


quadratique généralisé Ra. 


Théorème 2.6.1 

On note par DE, DE et DE les dérivées au sens de Fréchet en A,B et C' respectivement. 

Pour tout z E re(A, B,C), on a 

e L'opérateur Ra(z,:, B,C) : BL(B) — BL(B) est différentiable au sens de Fréchet et pour 
tout opérateur H € BL(B) sa dérivée au sens de Fréchet en A est donnée par la formule 


suivante : 
D1Ro(z, À, B,C)H = —-Ra(z, À, B,C)H Ra(z, À, B,C), 


e L'opérateur Ra(z, À,:, 0) : BL(B) — BL(B) est différentiable au sens de Fréchet et pour 
tout opérateur H € BL(B) sa dérivée au sens de Fréchet en B est donnée par la formule 


suivante : 
DB Ra(z, À, B,C)H — —2Ro(z, À, B,C)H Ra(z, À, B,C), 


e L'opérateur Ra(z, A, B,:) : BL(B) — BL(B) est différentiable au sens de Fréchet et pour 
tout opérateur H € BL(B) sa dérivée au sens de Fréchet en C' est donnée par la formule 


suivante : 


Do Ra(z, 4, B,OH — — Ra(z, À, B,C)H Ra(z, À, B,0), 


Preuve. 

Soient z € re(A, B,C') et À, B, et C trois opérateurs dans BL(B). On va montrer qu’il existe 
un opérateur qui est la dérivée au sens de Fréchet de RQ(z,-, B,C) en A. 

On peut écrire 


Ro(z, A+H,B,C)-Ra(z, A, B,C) = (*(A+H)+2:B+C)1-(A+2B+C) 
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On prend (2?(A + H) + 2B + C)-! comme facteur commun à gauche et (22A + 2B + C)"! 
à droite. On obtient alors, 


Ra(z, À + H, B,C) = Ra(z, 4, B,C) 
((A+H)+2B+0C)1|24+2B+C-(2(A+H)+2B+0C) (2A+2:B+0C) 


((A+H)+2B+0C)1(-2H) (PA +2B +0) 1 = -2Ro(z, A+ H,B,C)H Ro(z, 4, B,C), 


Ra(z, A+H,B,C)-Ra(z, A, B,C) = Ra(z, A, B,C)-2Ra(z, A+H,B,C)H Ra(z, À, B,C\. 
Ceci donne 


Ra(z, À + 4,B,C) - Ra(z, 4,B,0C) 
=. [Ro(z, 4, B,C) — 2Ro(z, A+ H,B,C)H Ro(z, 4, B, C)] HRo(z, A, B,C) 
= —2Ro(z, À,B,C)H Ra(z, 4, B,C)+2“Ra(z, A+ H,B,C)H Ra(z, A, B,C)HRa(z, A, B,C). 


Il est clair que 
zRa(z, 4, B,C)H Ro(z, A+ H,B,C)H Ra(z, A, B,C) = o(||H|°) 


pour ||H| proche de 0. 
Ceci signifie que D1 Ra(z, À, B,C)H = —-2?RQ(z, À, B,C)H Ra(z, A, B,C). 


De manière similaire, on obtient les autres résultats. 


2.7 Conclusion 


Dans ce chapitre, nous avons étendu les concepts de la théorie spectrale du cas généra- 
lisé à notre cas, afin de définir le spectre quadratique généralisé associé à trois opérateurs 
linéaires bornés. De plus, nous avons également développé les outils théoriques de la méthode 
spectrale quadratique généralisée. Enfin, nous avons généralisé certaines propriétés classiques 
mentionnées dans [I]. Ces généralisations nous ont permis de construire un nouveau type de 


convergence dans le deuxième chapitre. 
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3.1 Introduction 


Notre principal intérêt réside dans la détermination des valeurs spectrales associées aux 
trois opérateurs linéaires bornés À, B et C', mentionnés dans le premier chapitre (section 
1). Comme la résolution directe des problèmes spectraux est assez compliquée, car le calcul 
des valeurs spectrales exactes est presque impossible dans la majorité des cas, nous sommes 


contraints d'établir des estimations approximatives pour ces problèmes spectraux. 


Généralement, ces approximations sont les résultats exacts de calculs du spectre quadra- 
tique généralisé associé aux trois opérateurs linéaires bornés À, B et C' qui sont des approxi- 


mations des opérateurs À, B et C' donnés. 


Souvent. les opérateurs À, B et C' sont choisis, pour être éléments des suites d'opérateurs 
linéaires bornés (4,),(B,) et (Ch), convergeant dans un sens approprié vers À,B et C 


respectivement. 


Pour les raisons évoqués ci-dessus, on se concentre sur la construction d’une méthode 


approximative pour le spectre quadratique généralisé. 


Soit (B, || -||) un espace de Banach. Soient À, B et C trois opérateurs linéaires bornés dans 
BL(B), l’espace des opérateurs linéaires bornés de B vers B. Soient {4,},cN; {Bn}ren €t 


{Cn}nen trois suites dans BL(B). 


Pour bien construire notre méthode et pour déterminer la relation entre le spectre qua- 
dratique généralisé approché sp(4,,, B,,C,) et le spectre quadratique généralisé sp(A4, B, C), 
nous supposons que les suites d'opérateurs {A}, ; {Bnjnen €t {Cnnen Convergent vers 
les opérateurs À, B et C’ respectivement, dans un sens donné qui sera précisé ultéreurement. 
Aussi, il faut suivre le schéma d’approximation associé au problème spectral quadratique, 
développé dans le premier chapitre, à savoir, 


Propriété U : 


Si Àn € SP(An, Ba, Cn) et Àn — À ,alors À € sp(A, B,C). 
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Propriété L : 


Si À e sp(A, B,C'), il existe une suite À, € sp(A,, Br, Cr) telle que À, — À. 


3.2 Approximation du spectre quadratique généralisé 
au sens de la propriété U 


Dans cette section, nous allons vérifier la propriété U, mentionnée ci-dessus, qui justifie 
l’absence de la pollution spectrale pour les opérateurs bornés. Aïnsi, nous allons prouver la 
convergence du spectre quadratique généralisé approché sp(A,, B;, C), vers le spectre qua- 
dratique généralisé exacte sp(A, B,C'), en utilisant la convergence des suites des opérateurs 
{An}nen > LBn}nen: {Cn}nen n norme ainsi que la convergence en collectivement compact, 
vers les opérateurs À, B, C' respectivement. 
De plus, nous allons procéder à l’extension du nouveau type de convergence ( 7—convergence généralisée) 


développé dans le premier chapitre au cas quadratique. 


Le théorème suivant nous permet de démontrer notre résultat principal qui sera développé 


ultérieurement. 


Théorème 3.2.1 
Soient À, B,C! À, B et C dans BL(B), etzec. 
On définit les opérateurs Q(2) et Q(2) sur BL(B) par 


Q(2) = AB + C,, (3.1) 

Q(2) = À +2:B+ÛC. (3.2) 
Pour z un élément de re(A, B, C') tel que Q(z) et Q(z) vérifient la propriété 

ICO OILTEIES 3) 


on az ere(À, B,C). 


De plus, la norme de la résolvante quadratique généralisée de À, B et À, est donnée par 


(3.4) 
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et 
Ra(z) = [(Q()]7 = (24+2B+ 0) 


De même, si Q(z2) et Q(2) avec z € re(A, B,C) satisfont 
; 2 
ICO OL CINE (3.5) 
alors, on a z € re(À, B,C), en plus la norme de RG(z) dans BL(B) vérifie l'inégalité suivante : 


Rate) [1 + u -4@)Re(2)]"|] 


a C ao ah 


IRG(z)] < (3.6) 


Preuve. 

Pour montrer ce théorème, il suffit de prouver l'existence et la bornétude de l’opérateur 
RG (2), qui servent à montrer que z € re(À, B,C'. 

Soit z € re(A, B,C). Q(2) s'écrit 


Q(:) = Q(z) — Q(2) + Q(2). 
On prend Q(2) alors en facteur commun, ce qui donne 
QG) = Q(G)[1- (Q(2) -Q(2)) Ra()|, (3.7) 


avec RQ(z) = [Q(2)] 


— 
Aussi, pour établir que RG(z) = [Q(2)] existe, il suffit de montrer que les opérateurs 


[1 — (Q(:) —Q()) Ra(2)] et Q() 


sont inversibles. 
Il est clair que Q(z) est inversible et borné car z € re(A, B,C'), et on peut écrire RQ(z) = 
[Q(2)77". Maintenant, pour montrer que [I — (Q() — Q(2)) Ra(?)| est inversible et borné, 


en utilisant la condition (5.3) et d’après le théorème|2.5.2| on a 


—1 


[! - (QG) - 4) Ra] 


ce qui montre l’existence de l’opérateur inverse ci-dessus et que sa norme est bornée par 


1 
1 — | (Q(2) — Q(2)) Ra(2)| 


FÉACOMOLORE (3.8) 
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Cela implique que 


1 


Ra() = [62] = [1- (Q(2) - G(2)) Ro] Ra(). (3.9) 


De 3.8] la norme de R&(z) est majorée par 


I Ra(2)| 
. | 3.10 
ROIS Te -dG)RQ] L 


Donc, R&(z) existe et borné, ce qui montre que z € re(À, B,C). 


Ensuite, on montre que R&(z) est borné moyennant la condition B-5). 
Comme 
- (00 -00)RG) 
existe et est borné, d’après la série de Neumann, on peut écrire 


1 


1-(Q()-0())Rat] = > [(Q(2) — Q()) RaG)|” (3.11) 


ce qui nous permet d'écrire 


[ee] 


Rq() = Ro(9 Y [(Q(:) - 40) RG] 
On décompose la série en deux séries l’une des termes paires et l’autre des termes impaires. 
On obtient 
Ro) = Ro) > [(Q() - 40) Rata)]” 
+ > (Qt) - 4) Ra] Rae) 


Avec une simplification, on a 


k 


Ra(z) = Ra(z) [1 + (Q(:) - Q(:)) Ra(2)] 2 (QG - @()) Ra(:) | 


En appliquant la norme sur cette dernière égalité, on aboutit à 


I Ra) 11 < 1 Ra) [7 + (Q(2) — Q(2)) Ra(2)] | > | |((Q(2) - 00) Ra(c)) | Ps 


En utilisant le fait que la série 
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est une série géométrique de rayon 


(CROIS 


on conclut que cette série est convergente vers 


= | (Qt) 00) Rat) 


ce qui donne 


IR [1 + ][(@() - Ga) rat] 
1- [en -06) Ref 


Donc, les deux propriétés (3.4) et (3.6) sont vérifiées et le théorème est démontré. 


RG < 


Afin de prouver la convergence du spectre quadratique généralisé au sens de la pro- 
priété U, nous allons établir cette convergence sous les trois types de convergence à sa- 
voir, la convergence en norme, la convergence en collectivement compact et enfin la  — 


convergence généralisée . 


1, 4 nm 
3.2.1 Propriété U sous la convergence en norme (notée —) 
Le théorème suivant montre que la propriété U est vérifiée sous la convergence en norme. 


Théorème 3.2.2 
Soient {An }nen, {Bn}nen €t {Cn}nen trois suites d'opérateurs dans BL(B), vérifiant l’hypo- 
thèse : 


(H1) An + À, Br + B, CG +C, 


alors la propriété U (pour Ày € Sp(An, Ba, Cn), Vn € N telle que À, converge vers À dans C 
alors À € sp(A, B,C')) est satisfaite. 
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Preuve. 


La démonstration de ce théorème est basée sur le raisonnement par l’absurde. 
On suppose que À € re(A, B, C'). Comme l’ensemble résolvant quadratique généralisé re(A, B, C') 


est un ouvert dans C, il existe un scalaire à > 0 tel que 
QA={z2EeC:z2- A] <6} Cre(A,B,0\. 


De plus, la résolvante quadratique généralisé RQ(-) est analytique sur l’ensemble re(A, B,C) 


(voir le théorème 2.5.3), et on peut définir le réel 5 par 
B = sup{||Ro(2)]|] : z € A} < co. 


On cherche à montrer que l’ensemble Q est inclus dans l’ensemble re(A,, Bs, Cr). Pour cela, 
il suffit de prouver que pour z € (. on a z € re(A,, B,,C), ce qui revient à montrer que 


l'opérateur Ro, (2) = (274, + zB, + C,) ! existe et est borné. 


D'une part, en utilisant l'hypothèse (#1), on suppose, de plus que pour z € Q ils existent 


trois entiers naturels n1,n2 et n3 vérifiant 


1 


A A en ps 3.12 

PE (3.13) 
HSE} | 

IC - C, | < L Lu (3.14) 


On montre alors que 
IQ(:) — Q: (2) Ra) < TL 
Pour tout n > max(n1,n2,n3), on a la relation suivante : 
IQ) - QG RQ) = |[22(4 — 4,) + 2(B — B,)+(C - C,)] Ra(2)| 


Comme |Ra(z)|| < B et d’après (3.12) — (3.14), on a la majoration suivante : 


IQ) -Q IR < (4/14 - 41 + 1418 — Ball +11C - Cl) 8 


_ ST 
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Donc, d’après la condition (|[Q(z2) — Q,(z2)] Ra(z)|| < 1) l'opérateur 
Ro (2) = (27 An + 2Bn + Cr) 


existe et est borné. Par conséquent z € re(A;, By, Ch), et donc,  C re(A,, Bn, Cn) pour tout 


n assez grand. 


Comme la suite des valeurs spectrales approchées À, converge vers la valeur exacte À, et 
Â]n EN Cre(Ah, Br, Cn) pour tout n, À, n’est pas dans le spectre sp(A,, Bn, Ca), ce qui est 
en contradiction avec l'hypothèse de la propriété U (À, € sp(4,, B,,C,) pour tout n). Par 


conséquent, À doit être dans sp(4A, B, C). 


3.2.2 Propriété U sous la convergence collectivement compact (no- 
tée +) 
Dans cette section, on va prouver que la propriété U est réalisable sous la convergence des 
suites d'opérateurs À,,B,, et C, en collectivement compact. 
Dans la suite, on a besoin du lemme technique suivant pour démontrer notre théorème fon- 


damantal. 


Lemme 3.2.1 
Soient {N;},ens (Mn} 
dans BL(B). 


nen deux suites d'opérateurs dans BL(B), N,M et K trois opérateurs 
Si les suites N, et M, convergent vers N et M respectivement en collectivement compact, 
alors 


Preuve. 


Comme les suites {N,}.. et {M}. sont convergentes en collectivement compact, on 


nEN neN 
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Pour tout x € B, lim ICN — N,)x|8 = 0 et il existe un entier n, € N tel que l’ensemble 
Si = (J {Nix - Nx :x€B, [xs =1} 
n>ni1 


est relativement compact (définition de la convergence en collectivement compacte). 


Pour tout x € B, lim CM — M,)x|# = 0 et il existe un entier n2 € N tel que l’ensemble 


suivant 


S2= (J {Mix - Mr :x € B, [xls = 1} 


n>n2 


est relativement compact (définition de la convergence en collectivement compact). 


Comme À est un opérateur borné, nous savons que les ensembles X(S:) et K(S2) sont 


relativement compacts. Par conséquent 


CN — Nh)K(M — M) || — te CN — Nh)K(M — Mh)t||8 
L||B— 
— ; co IN — N,)K(Mzx — M,x)|s 
x||8= 


et pour y = K(Mx — M,x) € K(S2), on obtient la relation 


IN-NDAM- MI = sup IV Nvlls. 
yEK(S2) 


En passant à la limite, on a 


Bi JON N)AGM M) € sup Jim JON Nul 
yEK(S2) 


De plus, 
Tim [CV = N,)glls = 0. 


Finalement, on vérifie que 

Ti CN = N)K(M = M) = 0. 
De manière simililaire on obtient 

dim CM — M)K(N = N)I = 0. 


Cela montre le résultat (3.15) et donc, le lemme 3.2.1. 
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Dans le théorème ci-dessous, on vérifie la propriété U sous la convergence en collectivement 


compact. 


Théorème 3.2.3 
Soient {An}nen, {Bn}nen €t {Cn}nen trois suites d'opérateurs dans BL(B), 


où les opérateurs À,, B, et ©, vérifient l'hypothèse suivante : 
(H2) An + À, B, + B, Cr + C. 


Alors la propriété U {pour À, € sp(A,, B,, Ch), Vn € N telle que À, converge vers À dans C 
alors À € sp(A, B,C')) est satisfaite. 


Preuve. 
La démonstration de ce théorème est basée sur le même raisonnement du cas précédent. 


On suppose que À € re(A, B, C'). Comme l’ensemble résolvant quadratique généralisé re(A, B, C') 


est un ouvert dans C, il existe un scalaire à > 0 tel que 
Q={z2ec:l2-AÀA| <6}k Cre(A,B,0C\). 


De plus, la résolvante quadratique généralisé RQ(:-) est analytique sur l’ensemble re(A, B,C), 


et on peut définir le réel 5 par 
B = sup{||Ro(2)|] : z € A} < co. 


On cherche à montrer que l’ensemble Q est inclus dans l’ensemble re(A,, Bs, Cn). Pour cela, 
il suffit de prouver que pour z € on a z € re(A,,, B,,C,), ce qui revient à montrer que 


l'opérateur Ro, (2) = (224, + zBy + Ch) ! existe et est borné. 


D'une part, en utillisant l'hypothèse (7/2) ,on suppose, de plus que pour z € Q ils existent 


des entiers naturels n; avec à = 1,9 vérifiant 


[CA — 4) Ra(z)(A — A)]l < TEU) de. (3.16) 
I(B — Be) Ro(z)(A — À) - rs (317) 


I < BG+h "2 
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IAA) RQGNB BIS ET En (3.18) 
I(B—8,)Ra()(B — B,)|| < GENE n > m4, (3.19) 
I(A — An) Ra(z)(C — Ch) < ME D n > 3, (3.20) 
I(C — Cu) Ra(z)(A — 45) < BIG DEr n > n6, (3.21) 
(8 — B,) Ra(2)(C — CI < TENUE n > (3.22) 
IC — C)Ra(2)(B — B,)| < GET n > 8, (3.23 

NC G)RaGNC - CS ge 2e (3.24) 


En appliquant la norme, on obtient 


IQ) — QG) RGP] =" 1I(Q() — 0:62) Rata] (QG) — Q(2)) Re 
<_|I(@() — Qr(2)) Ra(z)(Q(2) — Qu (2) Ra (2). 


Pour tout n > max(n:,n2,n3, Na, 5, 6, 7, 8, 9), et d’après (3.16) — (3.24), on obtient la 


majoration suivante : 


IQ) - Q, (2) Ra(DP] < 8(I(C - C,) Ra(2)(C - C,)| 


+ |AI(A — 45) Ra(z A) + 1218 — Ba) Ro(z)(4 — 45) | 
B,)1 + 12 1(8 — B:)Ra(2)(B — B:)| 


(4 
(B 
(C— Ca) +12 I(C — Cr) Ra(z)(A — A) 


(C 
) 
+ 12°C — 4x) Ra(z) 
+ lI(A — 4x) Ra(z) 

( 


I(B — B,) RaG)(C — Ch) + AC — Cr) Ra()(B — B,)1) 


| 
L 
] 


ES 
TT Cr Er UT UT MT UÉTEÉT 


Donc d’après la condition 
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l’opérateur 


Ro:(z) = (A, + 2B, + Ci" 


existe et est borné. Par conséquent z € re(A,, B,, CA), donc Q C re(A,, B,,C,) pour tout n 


assez grand. 


Comme la suite des valeurs spectrales approchées À, converge vers la valeur exacte À, et 
]n EN Cre(A;, By, Ch) pour tout n, À, n’est pas dans le spectre sp(4,, B;, Ch), ce qui est 
en contradiction avec l'hypothèse de la propriété U (À, € sp(A,, B,, Ch)) pour tout n. Par 


conséquent, À doit être dans sp(4A, B, C). 


e » L L L L # # V, 
3.2.3 Propriété U sous la v — convergence généralisée (notée —) 


Dans cette section, on va introduire une nouvelle convergence, appelé v — convergence 
généralisée, où la propriété U sera toujours réalisable sous des conditions plus faibles que 
(H1) et (#2). Pour cela, on va faire une extension de la nouvelle notion de v — convergence 


généralisée. 


La définition suivante permet d’introduire la v — convergence généralisée. 


Définition 3.2.1 

Soient An, Bn, Cn, À, B,C des opérateurs dans BL(B). On dit que le triplet des suites d’opé- 
rateurs (An, Br, Cn) converge en v — convergence généralisée vers le triplet d'opérateurs 
(A, B,C'), et on note (An, Br, Cn) 5 (A, B,C') si les conditions suivantes sont satisfaites : 


(i) Les normes ||A,{|, [[B;]| and ||C;]| sont bornées indépendamment de n. 


(ii) || (An — À) Li] — 0, [ (Br — B)L,l] — 0, [ (Cr — ©) L,ll — 0, pour n tend vers l’infini 
où L, prendra toutes les suites AÀ,, B, et Cn,(Ln € {An, Bn, Cn}). 


(iii) || (An — À) L] — 0, | (Bx — B) L] — 0, | (Cr — ©) L] — 0, pour n tend vers l’infini 
où L prendra tous les opérateurs À, B et C, (LE {A,B,C}). 


Remarque 3.2.1 


Cette v — convergence généralisée sera développée pour le traitement du spectre quadratique 
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généralisé de la forme sp(A, B,I1 + D), où I est l’opérateur identité dans BL(B), et À, B,D 


sont trois opérateurs dans BL(B). 
On considère le problème spectral quadratique généralisé sous la nouvelle forme sp(A, B, 1+ 
D) tel que 
Trouver ÀEC et ue B—{0} tel que Au + ÀBu + (1 + Dju = 0. (3.25) 


Dans le théorème suivant, nous allons prouver la convergence du spectre quadratique géné- 


ralisé associé au probème (3.25), au sens de la propriété U, sous la  —convergence généralisée. 


Théorème 3.2.4 
Soient {An}nen, {Bn}nen et {Dunen trois suites d'opérateurs dans BL(B), 


où les opérateurs À,, B, et D, vérifient l'hypothèse suivante : 
(H3) (A, Ba, Dr) + (A,B, D). 


Alors la propriété U (pour An € Sp(An, Bn, 1 + D), Vn € N telle que À, converge vers À 
dans C alors À € sp(A, B,1I + D)) est satisfaite. 


Preuve. 
La démonstration de ce théorème est basée sur le raisonnement par l’absurde. 


On suppose que À € re(A, B,1 + D). Comme l’ensemble résolvant quadratique généralisé 


re(À, B, 1 + D) est un ouvert dans C, il existe un scalaire 0 > 0 tel que 
QA={z2EeC:|z2-A|<6}Cre(A,B,I+D). 


De plus, la résolvante quadratique généralisée RQ(:-) est analytique sur l’ensemble re(A, B,C), 


et on peut définir le réel & par 
a = sup{||Ro(z)|| : z € A} < oo. 


On cherche à montrer que l’ensemble Q est inclus dans l’ensemble re(A,, B;, Cr). Pour cela, 
il suffit de prouver que pour tout z € Q alors z € re(A,,, B,, 1 + D,), c’est à dire on montre 


que l’opérateur Ro, (z) = (274, + z2B, + C,) ! existe et est borné. 
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Pour tout z € (, on pose 
Qh() = A, + 2Bn + (1 + D) 


et 
Q(z) = A+z:B+(I1+D). 


D'une part, on va développer le terme [(Q(z) — Q,(z)) Ra(z)]”, en utilisant l’egalité pré- 


sentée dans le Théorèm 2.5.4(b), qui nous permet d'écrire RQ(z) sous la forme suivante : 
Ra(z) = 1 — 2 ARa(z) — zB Ra(z) — DRa(z). 
Ceci donne 


[(Q(z) — Qn(2)) Ra = (Q(z) — Qn(2)) Ro(z)(Q(z) — Qh(2)) Ro(z) (3.26) 
= (Q(2) — Qn(2)) [I - A Ra(2) — 2BRa(z) - DRo(z)| 
(Q(:) — Q,(2)) Ra() 
= {(Q(:) — Q,(2))Q(2) — (Q(2) — Qh(2))Qh (2) 
—  2(Q(2) - Q(2)ARa(z)(Q(2) — Qn(z)) — 2(Q(2) — Qn(2))B Ra(2) 
(Q(z) — Qn(2)) — (Q(2) — Qn(2))P Ra(z)(Q(2) — Qn(2)) Ra(2) 
D'autre part, du fait que ||A,||,||B,| et |1D, I] sont bornés, , il existe t1,t2,t3 > 0 tels que 


A — À, < &,11B — B,l| < t2 et | D — D,| < t3 pour tout n EN. 
D’après (#3), on a 


| (An = A) El, | CB = B) El | (D; a D) L, || eu 0, La € A Bh, Dh} 


et 


De plus, on suppose que pour z € Q, et to — max(f1,t2,t3), ils existent des entiers naturels 


n;,i = 1,18 vérifiant 


1 
So. 3.27 
IC JA] < 36a(r + [A+ (X + + Mol) 1 (3.27) 
He ne Fe 
n)An|| = 36a(r + JA 4)? N 7 N4, ; 
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IAA) x EN de 
IA ADE <> 

I(A— A) PI < 36a(r + |X2- x DES O TT) SL. 
(A ADD € 2e 

ICE — BA < 36a(r + | - : EME 
(BBA <> 

ICE — B)BI< 36a(r + [X - x DES) 
(B-BDE <> 

IE BDIS En: x ant 
(BBD S pp Meme 

(D — DA)AÏ < 36a(r + [X - x Ep 20 
(D DA < em 

ID DIAIE re x EE Na = us 
lD-D)BIS rpg 2m 

IDD) = Sn 


D| < 
LE gat++à+0a ” 
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(3.29) 


(3.30) 


(3.31) 


(3.32) 


(3.33) 


(3.34) 


(3.35) 


(3.36) 


(3.37) 


(3.38) 


(3.39) 


(3.40) 


(3.41) 


(3.42) 


(3.43) 
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1 
I(D-D)D< SE n> mue. (3.44) 


En utilisant les inégalités précédentes de manière similaire aux cas de preuves des théorèmes 


3.2.2 et 3.2.3, on montre que 


ICORIOL ONE" 


Pour tout n > max{n;, 1 < i < 18}, on obtient 


ICOBAOL ONESES 


Ceci montre que 


Ro,(z) = (24, +2B, +(1+D,)) 7! 


existe et est borné, et z € re(4,, B,, 1+D,). Donc, Q C re(4,, B,, 1+D,) pour tout n assez 
grand. 

Comme la suite des valeurs spectrales approchées À, converge vers la valeur exacte À et 
MX EN Cre(A,, B,,1+D,) pour tout n, À, n’est pas dans le spectre sp(A, B, 1 + D), ce qui 
est en contradiction avec l'hypothèse de la propriété U (X, € sp(A,, B,, I + D,) pour tout 
n). Par conséquent, À doit être dans sp(A, B, 1 + D). 


Remarque 3.2.2 


— Si le triplet (A,, B,, D,) converge en norme vers le triplet (A, B, D) 
(An; Bns Dn) + (A, B, D) 


alors le triplet (A4, B,, D,) converge en v— convergence généralisée vers le triplet (A, B, D) 
(An: Bn; Dn) + (A, B,D)). 


— Si(A», Br, Dn) + (A, B, D) et les opérateurs À, B et D sont compacts, alors 


(An; Bns Da) + (A, B, D). 
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3.3 Approximation du spectre quadratique généralisé 
au sens de la propriété L 


Dans cette section, on va prouver que la propriété L est réalisable, et on va établir un 
nouveau cadre théorique pour appliquer la méthode d’approximation spectrale généralisée 
au problème spectral quadratique. De plus, on va étendre la méthode d’approximation dé- 
veloppée dans le cas du spectre généralisé au spectre quadratique généralisé afin de prouver 
la convergence de ce dernier au sens de la propriété L, sous la convergence en norme et la 


convergence en collectivement compact. 


3.3.1 Transformation du problème spectral quadratique en pro- 
blème spectral généralisé dans l’espace B x B 


Dans cette partie, on va démontrer une équivalence entre la résolvante généralisée et la 


résolvante quadratique généralisée moyennant certaines transformations. 


Le problème du spectre quadratique généralisé se transforme grâce à cet outil en problème 


du spectre généralisé. 
Ainsi, 
— On note par 4 = B X B, l’espace de Banach associé à la norme suivante : 
VU = (ui,u) € #, [U|x = fulle + Mulls- (3.45) 


— On note par BL(Y), l’espace de Banach de tous les opérateurs linéaires bornés définis 


sur X dans lui même, et sa norme est donnée par 


VAEBL(x): [AÏ.= sup [AUl+. (3.46) 
IU Lx =1 


Le théorème suivant permet de montrer l’équivalence entre la résolvante généralisée et la 


résolvante quadratique généralisée. 


Théorème 3.3.1 
Soient À, B et C' trois opérateurs dans BL(B), et z € C\ {0}. 


Si l'opérateur C' commute avec Ra(zi) et Ra(z2), avec 4,22 € re(A, B,C') \ {0}, alors, 
Q(2) = A +2B +C est inversible si et seulement siz € re(T, S), 
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où, les opérateurs T et S sont dans BL(X), donnés par la représentation matricielle suivante : 


rs (7 2) s=( 6 a (3.47) 


Let O sont respectivement les opérateurs identité et nul dans BL(B). 


Preuve. 


Pour montrer ce théorème il faut établir que l’opérateur inverse 
QG) = (24+2B+0)" 
existe si et seulement si z € re(T,S) ,c-à-d, 
(T -28S) 


existe. 


On commence par montrer que si 
me 
Qu = (24+:B+0) 


existe, alors 
(T - 28)" 
existe. 
Tout d’abord, on montre que les opérateurs [Q(z2)] et (T — zS) sont bornés. 
Soit z € C\{0}. On a alors, 


ms (7154 … 


D'après l'écriture des opérateurs Q(2) et (T — 2S), on remarque que ces derniers sont bornés. 


En effet, par un calcul des normes on a 


2 
TI = max À ITills = max(|B1|3 + [Cle 1), (348) 
Si2 
et 
2 
(SI = IRIS 2 LSills = max(| As, 11e): (3.49) 


En utilisant la bornétude des opérateurs À, B,C', on prouve que l’opérateur (T — zS) est 


borné. 
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Maintenant, on suppose que Q(z2)7! 


existe et on va montrer que l’opérateur inverse de 
(T — zS) existe pour tout z # 0. 
D’après l'application des calculs d’inverse d’une matrice sur le block d'opérateurs (T — zS), 


on obtient 


zQ(z) T7  —-Q(2)"C 
| Q(2)7 na (3.50) 


En effet, il suffit de prouver que (3.50) est l’inverse à gauche et à droite de (T — zS). 
Pour l’existence de l’inverse à droite, en utilisant la multiplication des matrices, on a 


(T 28) 2Q(2)77  —-Q()7 2: PB 0 2Q(2)77  —-Q()7 
QG)" Q()'(B is —1 21 QG) Q( LB + zA 
| I —-(B+2:A)Q(2)7!C +CQ(2)7(B + zA) | 
O I 


D’après le résultat obtenu dans (a) du théorème 2.5.5, on obtient 


QT QC \_{r0 
FE | Q@7 Q47{(8+z:4) Jo 1} 
Pour l'inverse à gauche, on procède de manière similaire au traitement du cas précédent. En 
utilisant encore le résultat obtenu dans (a) du théorème 2.5.5, on obtient 


QG QC in (#1 00) ic B+:A C 

Q()1 Q()(B + 24) out a(8+:4) LI I 
_ 2Q(2)74(B + 24) + Q(2) 2Q()71C — 2Q(2)71C 

Q() (BE + 24) - Q(- \(B + :4 QC + Q(2) 7 (zB + A) 


oo 


Alors, l'inverse à gauche et à droite de l'opérateur (T — zS) existe si Q(z2) ! existe, et 


l'inverse de l'opérateur (T — zS) est donné par 


0) —Q()!C 
ms (02 on w rca ) ou 


Maintenant, on va montrer que si l’inverse de l’opérateur (T — 2S) existe alors Q(2) ! 
existe. En effet, pour tout z Z# 0, si (T — 2S) 7 existe, on peut l'écrire sous forme d’une 


matrice par blocs de la manière suivante : 


Ku Ki 
_:S | 3.52 
msn -(in ne) (3.59 
où X;j € BL(B), i,j = 1,2. On a alors 
RE TS BASA CL LE O 
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où 
| HU | pi 56) (3.54) 
(3.53) peut s’écrire sous forme du système suivant : 
(1) Ku(B+2A)-Ki2 =1I 
Une (3.55) 
(4) KanC + 2K3 0; 
Ainsi, à partir de l'équation (2) dans (3.55), on a 
Kia = —=KuC. 


En remplaçant X,2 dans la première équation du système (3.55) et en mettant -KX;; en 


facteur, on obtient 
(3.56) 


1 1 
=Ku(A + 2B + C) — zKuQ(z) — À, 


Ceci montre que l'inverse Q(z)"! existe. 


D'’apès les deux implications montrées ci dessus, on conclut que 


Q(2) = A +2B +C est inversible si et seulement si z € re(T,S), Vz Z 0. 


Le lemme suivant nous permet de montrer que si les suites d'opérateurs 4,, B, et C, dans 
BL(B) convergent en norme et en collectivement compact, alors les suites de matrices par 


blocs d'opérateurs T, et S, dans BL(Y), convergent aussi en norme et en collectivement 


compact. 


Lemme 3.3.1 
Soient T et S deux matrices par blocs d'opérateurs dans BL(X) défini comme dans (3.47). 


Soient S, et T, deux suites de matrices par blocs d'opérateurs dans BL(X) données par, 
Bree As 0) 
(395) s-(0 25) 


(i) Soit l'hypothèse (H1): A, + À, B, + B, Cy + C. Alors 


(H1): TST, 5, +S. 


60 


SomiA. KAMOUCHE UNIVERSITÉ 08 Mai 1945 GUELMA 


(äi) Soit l’hypothèse (H2): A, + À, B, + B, C,  C. Alors 
(H2): TL, ST, $, +S. 


.n CC Je ‘ 
Ici — et — désignent la convergence en norme et la convergence en collectivement compact 


respectivement. 


Preuve. 
(i) On va montrer que les suites T, et S, convergent en norme vers T et S respectivement. 
SoinenN.Ona 
Ta — T} = max(||B, — B| + ||Cx — C0), (3.57) 
et 
[Sn — SI] = max(|] A4 — 4], O1). (3.58) 


En utilisant (#1), on obtient [T, — T| — 0 et |S, — S] — 0. 
Cela prouve que les deux suites T, et S, convergent en norme vers T et S dans BL(X) 


respectivement, et (i) est vérifié. 


(ii) Pour montrer que T, et S, convergent vers T et S en collectivement compact BL(Y#) 


respectivement, on va prouver les points suivants : 


1. On montre que T, 2 T'et S, + S (4 désigne la convergence ponctuelle). 


Soient x = (41,2%2)E XetneN.Ona 
Tr — T)rfx = max(](Br — B)211l5 + (On — C}r2ll8, 1O Is), 


(Sn — S)xilx = max(|(A — A,)xills, OI), 


où © désigne l’opérateur nul. Selon l'hypothèse (#2), on a 


dira (Tr — Della = Jin, 1Sn — Sex = 0. 


n— 00 


Ceci prouve que T, et S, convergent ponctuellement vers T et S dans BL(Y#) 


respectivement. 
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2. Pour nr et n: dans N, on va prouver que les ensembles suivants 


Q = (] {Tor- Tr re X, xllx =1} 
n>no 
et 
NO = (Ù] {S:xr-Sx :re X, [xllx =1} 
n>n1 


sont des sous-ensembles relativement compacts de #. 
En effet, on peut écrire 


Q = (] {(Ba- Br +(Cn-C)rr : m1, € B, [xls <1, [xls < 1} x {O} 
n>no 

Q — U {(A» = A)x; :41 € B, xl < 1} X {0}. 
n>n1 


D’après (#2), les ensembles 


{(Ba — Bjr + (On — Cire : m1,70 € B, [ils < 1, Irells < 1} (3.59) 
et 


{An — Ajr : M1 € B, [tie <1} (3-60) 


sont relativement compacts dans B. 

Par ailleurs, il est clair que le singleton opérateur nul {O} est un ensemble re- 
lativement compact de B, et comme le produit de deux ensembles relativement 
compacts est un ensemble relativement compact, on obtient que (; et ( sont des 


ensembles relativement compacts dans #. 


D'après les points (1) et (2), on conclut que $, et T, convergent dans un sens collectivement 


compact vers $ et T respectivement. 


3.3.2 Propriété L pour le spectre quadratique généralisé 


On va montrer que les techniques utilisées dans le cas généralisé peuvent s'étendre aux 
matrices de blocs d'opérateurs, car les opérateurs T,,S,,T et S satisfont aux conditions 


utilisées dans l’étude du spectre généralisé. Maintenant, nous appliquons les techniques et 
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résultats du paragraphe précédent pour démontrer notre résultat principal à savoir que la 


propriété L est réalisable. Ceci sera prouvé dans le théorème suivant. 


Théorème 3.3.2 
Soient {An}nen, {Bn}nen €t {Cn}nen trois suites d'opérateurs dans BL(B). 
Soit À une valeur spectrale quadratique généralisée isolée dans sp(A, B,C). Supposons que 


l’une des hypothèses suivantes est vérifiée, à savoir 


(HA: 20 étA, SA BB GC 


où 


(H5): ÀZO, et A SA, B SB, Cn $C 


Alors, s'il existe une suite Àn € SP(A», Bn, Ch), pour n assez grand, alors À, converge vers À 


dans C. 
Preuve. 


Soient À € sp(A, B,C) et lun contour de Cauchy, qui sépare À de sp(A, B,C)\ {A}.Ce 


qui implique d’après le théorème 3.3.1 que À € sp(T,S). 


On définit la suite À, comme suit : 
Àn = int(l) Nsp(T,, Su). (3.61) 
On fixe r’ > 0 tel que la suite (A,)1en appartient à l’ensemble : 
E:={2ec: |1-2l<7r}. 


Maintenant, on va prouver que À, — À, dans C, avec À € sp(T,S). Pour cette raison, il suffit 
de montrer que toute sous-suite (À,(x) de la suite (An)nen, Converge vers À. 

Soit (A,(#)) une sous-suite qui converge vers À £ À. Puisque la propriété U est réalisable sous 
les convergences en norme et en collectivement compacte, on a À € sp(T,S). Mais, he fet 


sp(T,S)NE := {À}. Par conséquent, À = À, et donc À, — À, ce qu'il fallait démontrer. 
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Remarque 3.3.1 
Sous les conditions (HA) et (H5), on peut montrer que la propriété U est aussi réalisable, en 
suivant de manière similaire la procédure précédente qui nous a permis la réalisation de la 


propriété L. 


3.4 Conclusion 


Ce chapitre représente le résultat essentiel de notre recherche, à savoir une extention assez 
rigoureuse du cas généralisé au cas quadratique généralisé. Ceci a nécessité une généralisa- 
tion de certaines propriétés et techniques classiques, qui nous ont permis la construction et 


l’adoption d’un nouveau type de convergence appelé —convergence généralisée. 


Ce nouveau concept donne une approximation pour le spectre quadratique généralisé de 
(A; B;T + D), où les opérateurs sont seulement bornés. Ainsi, on a montré la réalisation de 
la propriété U sous les trois types de convergence ( la convergence en norme, la convergence 
en collectivement compact et la — convergence généralisé). De plus, en se basant sur le 
bloc d'opérateurs, notre problème du spectre quadratique généralisé, moyennant certaines 
techniques, est transformé en un problème du spectre généralisé dans un espace produit, où la 
propriété L s’avère réalisable sous la convergence en norme et la convergence en collectivement 


compact pour toute valeur spectrale non nulle. 
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4.1 Introduction 


Dans ce chapitre, on va appliquer la méthode du spectre quadratique généralisé au pencil 
quadratique de l’opérateur de Schrôdinger. Ensuite on réalise quelques tests numériques pour 


valider nos résultats théoriques, afin de confirmer l'efficacité de notre méthode. 


4,2 Construction d’un exemple de validation 


On va étudier "le pencil quadratique" de l’opérateur de Schrôdinger, défini dans L?([0, +oo|) 


par l’équation différentielle du second ordre suivante : 


(P) { ci + QAU(x) — y =0, x € [0,+oo, 


où, V,U sont des fonctions complexes et À un paramètre spectral. 

1 
Dans toute la suite, on choisit les fonctions V(x) = x? et U(x) — 3 de sorte que le problème 
(P) s'écrit : 


(P') { . no — + X)y = 0, 


Cependant l’étude approfondie du spectre quadratique généralisé, nécessite l’introduction de 


l'opérateur non borné $ : D(S) — L?([0,-+co[) donné par 
Sy = —-y" +x°y, x € [0,+ool, (4.1) 
où, 9 est l’oscillateur harmonique défini sur le domaine 
D(5) = H?([0,+o0[) n Le € L?([0, +oo|) : A 2 lp(x)?dr < +00} 


dans L?([0, +oof). 


Pour a €[0, +, on introduit l'opérateur $, donné par 


Sy = =" + y, y € D(S;) = H°(0, a), 


Remarque 4.2.1 
Pour calculer le spectre de l'opérateur S, on a besoin d'introduire la relation Va > 0,  sp(Ss) € 
sp(S) qui permet de déduire le spectre de l'opérateur S à partir du spectre des opérateurs S, 


définis sur des intervalles bornés de la forme (0, a]. 
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Par conséquent, d’après la remarque ci-dessus, on définit l’opérateur Laplacien noté L,, par 


L:D(Lp) — L?(0,a), 


4 


y + Lyy = -Y, 


où le domaine D(L,) de L, vérifie 


D(Ly) = {y € H°(0,a), y(0) = y(a) = 0} = H5(0, a). 


Le théorème suivant permet de montrer que l’opérateur Laplacien L, est inversible. 


Théorème 4.2.1 
L'opérateur L, est un opérateur inversible et son inverse est l’opérateur noté À défini sur 


L?(0, a) par 
Ay(z) — ROTON (4.2) 


où K, est le noyau de Green donné par 


Un EE Pr 


K(x,t) = a 43 
me ER Se 


Preuve. Voir [5], [10]. 


En utilisant le théorème précédent, on peut transformer le problème non borné (P/), dé- 
fini ci-dessus, en un problème spectral quadratique associé à trois opérateurs bornés, en 
appliquant l’opérateur inverse du Laplacien. Cette transformation nous permet de définir les 
opérateurs bornés B et C, qui agissent tous les deux sur l’espace L?(0, a), selon les deux 


formules suivantes : 
By(e) = — | Ka(x,Dy(bdt, (4.4) 
0 
et 


Cy(x) — te) + f° Ka(x, DE y(Oat (4.5) 
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Le théorème suivant montre l’équivalence entre le problème des valeurs propres (P/) qui 
est donné par, 


Trouver u = À— ECet y € D(S,) — {0} tel que 


et le problème de valeurs propres quadratiques généralisées associées aux opérateurs À, B et 


C. 


Théorème 4.2.2 
Le problème de valeurs propres (P!) est équivalent au problème de valeurs propres quadratiques 


associé au problème suivant : 
(Q) : Q(Ü)y = X'Ay + ÀBy + Cy = 0. (4.6) 


Ainsi, À est une valeur propre quadratique de S, si et seulement si À est une valeur propre 


quadratique généralisée de trois opérateurs bornés (A, B,C). 


Preuve. 


Soient À une valeur propre quadratique de 5, et y € D($,) le vecteur propre associé à À. 
Par conséquent 


Say = (À - dy. (4.7) 


En explicitant S,, on peut écrire l'équation 7 sous la forme suivante : 
Loy + 2? y — ÀAy + Xy = 0, (4.8) 


où L, l'opérateur Laplacien. 


En appliquant l'opérateur inverse À de l’opérateur L, à l'équation on obtient 
XAy + By + Cy = 0, (4.9) 
où 


Ay(x) = | "Ke, Du(é)dt. (4.10) 


68 


SomiA. KAMOUCHE UNIVERSITÉ 08 Maï 1945 GUELMA 


By) = = RAOTON (4.11) 


Culz) = y(z)+ [ke yba, (4.19) 


Par conséquent, À est une valeur propre quadratique généralisée de (A, B,C) et y est le vec- 


teur propre quadratique généralisé associé à À. 


Inversement, soient À une valeur propre quadratique généralisée de (A, B,C), et y € 


L?(0, a) \ {0} le vecteur propre associé à À. On a alors 


X°Ay + XBy + Cy = 0. (4.13) 
Par conséquent, 
XAy— X | | Kalx, Oy(tdt + y(x) + Î Ke, O6y(6)dt = 0. (4.14) 
Ceci est équivalent à 
X'Ay — XAy + y + A(x°y) = 0. (4.15) 


Comme À est un opérateur linéaire, on a 
y = AOy — Xy— x y), (4.16) 


pour y € D(L,), et donc, (Ay — y — x?y) € L?(0, a). 
En appliquant L, à l'équation on obtient 


L,y = y — Xy — x°y, (4.17) 


ce qui implique que 


ay = Ày — Xy — L,y. (4.18) 
De plus, on a 
L'abtodr < Plus + Millie + Matos < +00 
On conclut que y € D(S,) \ {0}, et qu’on a 
X'Ay + ÀBy + Cy = 0. (4.19) 
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Par conséquent 
L,(WAy) + Lo(ABy) + Ly(Cy) = 0, 
ce qui s'écrit 
Say = (À — X)y. 
Ainsi, À est une valeur propre quadratique de $,. 


Ce qui achève la preuve du théorème. 


4.3 Traitement numérique 


Dans cette section, pour établir des tests numériques nous allons étudier un problème 
spectral dont les valeurs propres sont connues. Nous signalons que ce problème a déjà été 


étudié dans les travaux et [IUI. 


Soit le problème des valeurs propres de l’oscillateur harmonique S' définie par (4.1). 


Ainsi, pour u € C et y € D(S) \ {0} on a 


Sy = y, (4.20) 


qui s'écrit sous la forme suivante : 


ns [ —y"+2y = uy, x € [0,+oof 
(P À y(0) = 0. 


On peut remarquer que le problème (P”) est équivalent au problème (P’) mentionné dans 
la deuxième partie de ce chapitre, ce qui implique que les valeurs propres de l’oscillateur 
harmonique $, peuvent être exprimées en fonction de À par l'équation u = (1 — X?). Ces 
valeurs propres sont bien connues et sont données par l’ensemble {An — 1, n € N*}. 


En particulier, si y prend la valeur 4n — 1 alors on peut déduire de l’équation y = (À — X?) 


VI6n — 5 


9 ? 


que les valeurs propres associées au problème (P’), sont données par À — : Et 
pour n € N:. 


D'après les résultats obtenus dans [I], on peut définir le spectre sp(S) par 


sp(S) = UÙ sp(Su), (4.21) 


a>0 
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où, l'opérateur S, est défini par 
Sa = —ÿ" +2°y, y(0) = 0, y(a) = 0. (4:22) 


Pour cela, il suffit d'étudier le spectre de $,. 


Soit le problème des valeurs propres S,y = uy , ce problème est équivalent à 
Trouver u € C et y € D(S,) \ {0} vérifiant 
"+ ay = py = (À jy, 
y(0) = ya) = 0. 
Dans la suite, nous allons utiliser trois techniques d’approximation, à savoir la méthode 
de Nystrôm, la méthode de Sloan et la méthode de Kantorovich, pour construire des suites 
d'opérateurs {A},en-: {Bn}new- + {Cn}nen- qui seront des approximations des opérateurs 


A,BetC.! 


On note par Æ, la subdivision de l'intervalle [0, a] donnée par 
En = ha =, 2; =(5—1)hn 1<j<nn2>2} (4.23) 
nm — 


4.8.1 Approximation par la méthode de Nystrôm 


La méthode de Nystrôm pour estimer les opérateurs intégraux À, B et C'est utilisée à cause 
de sa convergence en collectivement compact, comme établi dans la référence [2]. Grâce à 
cette approximation, nous pourrons définir les suites d'opérateurs {Ah}, ; {Bnhnenw+ €t 


{Cn}nen- de la manière suivante : 


Any(x) = > K,(t,ti)y(ti), (4.24) 
Euls) = — É K,(x,ti)y(ti), (4.25) 
Chute) — ) + + SK wK,(x,t:)Ey(t), (4.26) 


où, 


les poids w; vérifient, w1 = wy:1 = 2 et w,; = h pour 1 <i<n. 
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Par conséquent, lorsque on remplace x par t;, pour 1 < j < n +1, on obtient le problème 


spectral matriciel quadratique, issu de la discrétisation de 4,B et C, 


X M av + XnMpv + Mov = 0, (4.27) 
où, les matrices 
MA(i, 5) nn Kat; te); (4.28) 
Mi, j) —= —WiKa(t;,ti), (4.29) 
Mo(i,j) = n(é,5) +uiko(t; ti), (4.30) 


1h11 est la matrice d'identité d’ordre n + 1 et v est le vecteur propre quadratique. 
Dans la suite, on va utiliser la fonction "polyeig" dans MATLAB, pour calculer les valeurs 


quadratiques associées aux matrices (M1, M3, Mo). 


4.3.2 Approximation par la méthode de Sloan 


Dans cette section, on utilise la méthode de Sloan, pour construire les suites d'opérateurs 
{An hnen {Bn}new € {Cn}new. De plus, on montre que ces suites vérifiént les conditions 


de la v—convergence généralisée, mentionnée dans le deuxième chapitre. 


Tout d’abord, on rappelle que l’opérateur C, définie dans la formule (4.12), s'écrit comme 
suit : 


CDR T 


où 1 est l’opérateur identité de BL(B) et 
Dy(x) = . K, (x, Ey(b)dt. (4.31) 
0 


Ensuite, nous appliquons la méthode de projection de Sloan sur les opérateurs À, B, et D, 


pour construire leurs approximations définies par 
AS = Ans, B5 = Bn,, D° = Dr, (4.32) 


où 7, est l’opérateur de projection. 
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Le théorème suivant montre que les suites d'opérateurs obtenues par la méthode de pro- 


jection de Sloan, vérifient les hypothèses supposées dans nos études théoriques (chapitre 02). 


Théorème 4.3.1 
Soit n € N. Les suites d'opérateurs A$, B$ et DS convergent en v— convergence généralisée 


vers les opérateurs À, B et D respectivement. 


Preuve. 


Au début, on va montrer que les suites A$, B$, D$ sont bornées pour tout n € N° 


On note par LS € 4 B?: Der et la norme de LS vérifie 
PAST AE 


Comme l'opérateur L est l’un de trois opérateurs À, B, D, (L € {A, B, D}) et l'opérateur de 


projection 7, sont bornés, on conclut que LS est borné. 


Ensuite, on va montrer que les normes || (AS - A) L||, || (B$ — B) L|| et 
] Ge — D) L]|| tendent vers 0 quand n tend vers l'infini. 
Soit n € N°. On a 
I (A$ — A) LI = |A (mx — 2) LI 
Comme L est un opérateur intégral avec un noyau régulier, il est compact. De plus ñ,x — x, 
Vx € L?(0, a), et donc, cette norme || (AS — A) L]|| tend vers 0 quand n tend vers l'infini. 


De manière similaire, on montre que || (B$ — B) L|| et || (D$ — D) L]|| convergent vers 0. 


Enfin, on va montrer que les normes || (AS — A) LS|| ; || (B$ — B) LS] et || (DS — D) L>|| 
tendent vers 0 quand n tend vers l'infini. 


Soit L$ € {AS, B° pal En reconduisant les arguments précédents, on a alors 
I CAS — A) LÊI = IA (a — D) Lrall, 


CBS — B) LÊI = 18 (mr, — 7) Lrall, 


I (DS — D) LËI = D (ma — 1) Lral, 
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qui convergent vers (0. 


En se référant à ce qui a été établi ci-dessus, A$,B$ et DS convergent alors en 7 — 


convergence généralisée vers À, B et D respectivement. 


Maintenant pour obtenir un système matriciel, on va procéder à la discrétisation des opé- 
rateurs À,B et D. 
A l’aide de la même subdivision définie par (4.23) et en appliquant la projection de Sloan, 
pour transformer le problème de dimension infinie en un problème matriciel (voir [2]), on 


obtient 


où, 


r)= | Kia de(bdt, 8j) = | Kat, De(Dét, (0 = | Ka De; (Da, 
(4.33) 
{e; FE sont les fonctions chapeaux définies par les formules suivantes : 
pou 2<j<n—l,ona 


it 
h 
GA) PAT DÉTETs, 


si ji < t < t;, 


0 par ailleurs 


S1 
ex(t) h de = > 12; 
0 par ailleurs 


L— tn : 
ab=l 5 StriStSt, 
0 par ailleurs 


Le problème spectral matriciel quadratique, issu de la discrétisation de À, B et D en utilisant 


la méthode de Sloan s'écrit 


X Ms Un + AnMpsün + Mpsün + lun = 0, (4.34) 
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où À, est la valeur propre approximative, u, le vecteur propre associé à À, et les matrices 


Mis, Ms et Mhs sont données par 


Mas (i, 5) — @j(ti), (4.35) 
Mps(i,j) = Bj(ti), (4.36) 
Mo) = ft) (1.37) 


On rappelle que (4.35 — 4.37) seront associés à la fonction "polyeig" dans MATLAB, pour 


calculer les valeurs propres associées au triplet (Mps + 1, Mzs, M4s). 


4.3.3 Approximation par la méthode de Kantorovich 


Avec la même subdivision (4.23), on va construire les suites d'opérateurs {A4}, {Bn}nen 
et {D;},enw. Ensuite, on va appliquer la méthode de la projection de Kantorovich aux opé- 


rateurs À,B et D donnés par les mêmes formules présentées dans la méthode de Sloan à 


j | Klx, by(tdt. (4.38) 
= — Î "Klx, ty(t)dt, (4.39) 
Î | Ka(x, OEy(6)dt. (4.40) 


Rappelons que la méthode de Kantorovich est caractérisée par sa convergence en norme 


(voir [I]). 


Commençons par appliquer la projection de Kantorovich sur les opérateurs À, B, D, pour 


définir les opérateurs approximatifs AK, BK, DK (voir [2]). On obtient 


AKz(t) = D Arte => [Ka KA Je(s)ds| e,(t) alé > aje;(t 
sa En Em os 


où &;, 5; et y; sont donnés par (4.33). 
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Ainsi, le problème spectral quadratique sera réécrit sous la forme d’un problème de valeurs 


propres quadratiques de dimension finie : 
XD ae(t) + AD Bje;(t) + D ye;s(t) + x(t) = 0. (441) 
j=1 j=1 j=1 


En multipliant (4.41) respectivement, par X,(t;,s), —K,(t;,s) et K,(t;, s)s?, puis en intégrant 


sur [0, a], on obtient le système suivant : 


VI Pi, ja + À Pi(i, 5)85 + D Pa(i, j)y; + @i = 0, 


j=1 j=1 j=1 

X D_P(i,j)a; ca ÀD Pi, j)B; L = D Pi, j)y + Bi = 0, 
j=1 j=1 j=1 

VS P(i,j)a; + 1 Ps, 5)B; + S_ Pali, 5); + + 0, 
j=1 j=1 j=1 


où, les matrices P, P, et P3 sont données par 
Pi(i,3) = [| Kalti,s)e;(s)ds, 
0 
PA(i,5) = | -Kalts)e;(s)ds, 
0 


Pt) a Kit, s)s°e;(s)ds. 


Cela nous permet d’obtenir le système sous la forme matricielle par blocs suivante : 


P, 0 0 (02 0 P, 0 (02 Th 0 P, (02 0 

A2 | Px OÙ 0 HA TE À DILE |A) 0 & à B|=|0 

P3 0 0 e 0 P3 0 e DM 0D'ASE e 0 
(4.42) 


Le système (4.42) sera résolu en utilisant la fonction "polyeig" dans MATLAB, ce qui nous 


permet de calculer les valeurs propres quadratiques. 


4.4 Tests numériques 


Dans notre traitement numérique, on introduit l’argument suivant : 
An est dite € — valeur propre acceptable, si d(À,,sp(S)) < €, € > 0. (4.43) 


On note par d la distance entre À, et le spectre de l’opérateur oscilation harmonique $,,. 
Cette condition nous permet d'obtenir la propriété U numériquement et ainsi, le problème 
de la pollution spectrale est résolu. 


On note par 
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— VPp(n), le nombre de valeurs propres approximatives obtenues pour n > 1. 


— VP«(n), le nombre de £ — valeurs propres acceptables obtenues pour n > 1. 


= Ph = |VPan)/ Vlr] Yo 


Pour illustrer les résultats numériques obtenus par les trois approximations (Nystrôm, 
Sloan et Kantorovich), on a procédé aux tests numériques pour différentes valeurs de n avec 


e et a fixés. 


Les résultats obtenus par la méthode de Nystrôm sont présentés dans les tableaux. 


Tableau 4.1 concerne le cas : € = 107?, a = 1. 


mn VP» VPae PA 
100 198 166 83% 
200 398 374 93% 
300 598 570 95% 
400 798 764 95% 
500 998 974 97% 


TABLE 4.1 - Pourcentage des valeurs propres € = 107? — acceptable d’après [1.43] 


ré 
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Pour € = 107* et a = 1, les résultats numériques sont présentés dans le tableau 4.2. Pour 


mn VPp VPae PA 
100 198 20 10% 
200 398 106 26% 
300 598 232 38% 
400 T98 400 50% 
500 998 632 63% 
600 1198 860 71% 
700 1398 1074 76% 
800 1598 1276 79% 
900 1798 1466 81% 
1000 1998 1682 84% 
1200 2398 2074 86% 
1400 2798 2458 87% 
2000 3998 3678 92% 
2500 4998 4688  93.8% 


TABLE 4.2 — Pourcentage des valeurs propres € = 107* — acceptable d’après [43] 


les résultats obtenus par la méthode de Sloan avec a = 1 et £ = 1072, on a le tableau 4.8. 


nm VPm VPae PA 
100 195 165 84.62% 
200 395 365 92.41% 
300 567 567 95.29% 
400 795 763 95.97% 
500 995 967 97.19% 


TABLE 4.3 — Pourcentage des valeurs propres € = 107? — acceptable d’après [4.43] 
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Tableau 4.4 représente les résultats obtenus par la méthode de Kantorovich pour a = 1 


et € — 10 ? 
n  VPw  VPe PA 
100 196 166 8.69% 
200 396 366 92.42% 
300 596 068 95.30% 
400 796 764 95.98% 
500 996 968 97.19% 


TABLE 4.4 - Pourcentage des valeurs propres € = 107? — acceptable d’après [1.43] 


Nous comparons notre méthode avec la méthode des différences finies (FD) pour mieux 
montrer son efficacité et exposer le problème de pollution pour (FD). Nous utilisons la même 
grille Æ, définie ci-dessus. On note u; = u(x;) les valeurs approchées de u dans le noeud 
x;. En utilisant le schéma centré en trois points de la dérivée seconde, on obtient le système 


suivant 


Nous avons fait des tests numériques en utilisant la méthode des différences finies (FD) 
et l’approximation quadratique généralisée du spectre (GQS) avec la méthode de Sloan. On 
prend différentes valeurs de n et on fixe les valeurs a = 1 et £ = 10 ?. 

Tableau 4.5 propose de comparer les deux méthodes (FD) et (GQS) qui sont utilisées dans 
le calcul des valeurs propres quadratiques approchées. 


TABLE 4.5 - Pourcentage de la valeurs propres 107? — acceptable. 
n  nb-app-val nb-acc-val-FDM  nb-acc-valkGQS  Perc-FD%X  Perc-GQS% 


100 195 0 165 0% 84.62% 
200 395 6 365 1.5% 92.41% 
300 967 6 967 1% 95.29% 
400 795 6 763 0.75% 95.97% 
500 995 8 967 0.80% 97.19% 


4.5 Conclusion 


D'après les résultats obtenus dans nos tests numériques, on remarque que si la valeur de n 


croit, le nombre de valeurs propres acceptables augmente également, ce qui signifie que ces 
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valeurs convergent vers les valeurs propres exactes. À cet égard, l’application de la méthode 
d’approximation du spectre quadratique généralisée avec les trois méthodes (Nystrôm, Sloan, 
Kantorovich) permet d’éviter le problème de la pollution spectrale, apparu dans l’approxi- 
mation spectrale des opérateurs non bornés. 

Nos tests numériques ont montré l'efficacité de la méthode d’approximation du spectre qua- 
dratique généralisé basée sur les méthode de Nystrôm, de Sloan et de Kantorovich par rapport 
aux méthodes classiques. L'augmentation en pourcentage des valeurs propres de £—-acceptable 
dans notre méthode d’approximation du spectre quadratique généralisé permet d'éviter la 


pollution spectrale. 


Les résultats obtenus dans le dernier tableau montrent que lorsque la valeur de n est 
suffisamment grande, le pourcentage de valeurs spectrales acceptables augmente, et que les 
résultats obtenus par notre méthode sont bien meilleurs que ceux donnés par la méthode 
des différences finies. Ces résultats montrent que la méthode d’approximation spectrale évite 
la pollution spectrale, contrairement à la méthode des différences finies qui constitue une 
pollution spectrale lorsqu'elle est utilisée pour approximer des problèmes aux valeurs propres 


associées aux opérateurs non bornés. 
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Conclusion Générale 


La recherche présentée dans cette thèse a été consacrée à l’approximation du spectre qua- 
dratique généralisé des opérateurs bornés. Nous avons réussi à établir la propriété U et la 
propriété L, démontrant, ainsi, que le problème de la pollution spectrale n'intervient pas sur 
l’approximation des problèmes spectraux quadratiques associés aux opérateurs bornés. 

Ce résultat significatif nous permet de résoudre le problème de la pollution spectrale 
qui peut survenir lors de l’approximation spectrale d’un opérateur non borné en utilisant la 
méthode du spectre quadratique généralisé. 

Au cours de cette recherche, nous avons ouvert la voie à plusieurs perspectives futures. 
L'une de ces perspectives consistera à développer une méthode itérative visant à déterminer 
la meilleure valeur propre quadratique (maximale ou minimale), ce qui revêt une grande 
importance pour d’autres domaines de recherche. 

Une autre perspective porte sur la construction de l’opérateur de projection associé à la 
résolvante quadratique généralisée. Enfin, nous chercherons à développer une notion pour le 
pseudospectre quadratique généralisé. 

Cette thèse représente une avancée significative dans notre compréhension de l’approxima- 
tion du spectre quadratique généralisé et ouvre la porte à de futures recherches passionnantes 


dans ce domaine. 
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